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Fundamental Aspects Flow 


Part The Sphere Wake Steady Laminar 


review made the experimental evidence 
and theoretical considerations evolving from studies 
the wake formed aerodynamically-smooth 
sphere moving steady rectilinear motion through 
turbulence-free stream. The roles complex fluid 
dynamic phenomena which boundary layer 
separatioa, vorticity transfer and those characterized 
the Strouhal Number, well less evident fine 
processes are presented, with reference being 
made pertinent two-dimensional cylinder studies. 


URRENT studies the various aspects solids-gas flow 
are characterized growing attention the 
fundamental approach—a trend which many respects 
pical the changes which are occurring broader 
front all phases chemical engineering. One one, 
older concepts the science are being re-ex xamined 
the light the recent experimental evidence. Undeniably 
useful these concepts have been the past some 
them obviously require complete reappraisal, they 
longer appear account adequately for the increasing 
complexities modern technology. Similarly, the con- 
tinued usefulness friction factors, mass transfer co- 
and countless other semi- proportion- 
alities may many cases questioned, since they cer- 
tainly fail prov ide proper appreciation the under- 
mechanisms. extreme cases, particularly when 
coincident with ill use the dimensional 
analysis, such approach has been known result 
serious misconceptions. 


Increasing use the principles fluid mechanics 
the interpretation chemical engineering 
doing much improve the the basic 
phenomena. The necessity introducing large numbers 
assumptions for the solution complex 
problems being minimized the 
ability computing devices, while new mathematical 
methods approximation provide least tentative solu- 
tions even the most intractable non-linear 
equations. Finally, better instrumentation and the develop- 
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ment powerful experimental methods are contributing 
greatly the refinement the physical measurements. 

Nowhere are the beneficial effects these various 
improvements felt more strongly than the field 
solids-gas flow. There remains, however, the necessity 
re-assessing the old experimental evidence the light 
the recent findings. The function this series articles 
present critical analysis the more significant 
studies concerned with the various individual phenomena 
which, total, constitute the multiparticle solids-gas sys- 
tem. Starting with idealized single-particle behavior, each 
additional complexity will considered turn, laying 
the groundwork for subsequent discussion the com- 
plete problem. Part introductory concepts 
fluid mechanics were reviewed and the aerodynamics 
single smooth sphere moving very slow steady 
motion turbulence-free stream were des- 
the particle velocity relative the fluid 
increased, the momentum transfer becomes dependent 
the very complex phenomena the region wake, and 
these are discussed the present article. 


PROPERTIES THE WAKE 


Sphere Wake Configuration and Shedding Frequency 


The observation the wake sphere fraught 
with difficulty, and very often phenomena seemingly 
clearly established one observer will contradicted 
interpret even the stereoscopic photographs this region, 
and the results appear susceptible such factors 
the geometry the equipment and the turbulence the 
stream. 

illustration the conflicting evidence afforded 
the reports the first appearance the vortex ring 
which ushers the wake regime, and which probably 
comes into being the sphere boundary layer the rear 
stagnation the most downstream point the 
body the vicinity Reynolds Number ten. Care- 
ful observations this region convinced Williams that 
this did not occur even Reynolds Number 720, but 
Nisi and Porter clearly showed this occur 
8.14 photographing the smoke-filled flow air about 
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Figure 1—Sphere wake dimensions function the 
Reynolds Number (Taneda). 


supported sphere illuminated with powerful plane 
beam light. similar but possibly more refined tech- 
nique adapted water system was recently employed 
Taneda which showed the initial wake formation 
occur 24. Great care was taken eliminate the 
possible effect the wall, and Taneda felt that the sphere 
support rod, which can never done away with 
studies this type, exerted only minor influence 
the wake properties. The excellence the photographs 
obtained allows confident measurement the dimensions 
the wake and its development for the region 
300. and the results, which are for multi- 
particle studies, are given Figure They could not 
compared with the measurements and Porter, since 
the have failed mention the sphere size corre- 
sponding the plots eddy position which they had 
given. 

the Reynolds Number increases, the fixed ring 
vortex grows size and decreases stability and 
Reynolds value about 130 its 
begins oscillate. The oscillations become more severe 
Reynolds value referred the lower critical 
Reynolds Number detaches from the 
adhering the body the sphere and streams away 
velocity which less than the fluid relative velocity. 


Some that this occurs Reynolds Number 
200 and others 1000 but more likely first occurs 
value 500 for smooth sphere 
moving steady motion through turbulence-free and 
essentially infinite fluid. new ring forms the unstable 
one detaches, that there periodicity injected into 
the wake flow this Reynolds Number region (7), 
remarks that there should periodic increases 
and decreases the drag coefficient corresponding the 
alternate growth and away the rings, thus 
explaining the fluctuating fall velocities sphere noted 


There unanimity the fate the ring follow- 
ing detachment, but probably changes the flow 
observed that periodic shedding 
vortex rings existed only within the Reynolds Number 
range 500 1000, and Winny recorded photo- 
graphically the formation intertwining helical vortex 
filaments the range 2000 Ermish Rey- 
nolds Number 8400 reported the formation two 
vortex trailing the rear the sphere forming 


3000 4000 6000 10,000 


Figure 2—Sphere Strouhal numbers function the 
Reynolds Number 


interlaced spirals. Von Karman noted that ap- 
proximately 10,000 there existed two vortex filaments 
which were displaced 180° from each other, and whose 
position generation alternated irregularly. Lunnon 
points out that two such spirial threads indeed exist, 
they must act the sphere cause move 
slightly curved and spiral path. 

Perhaps the best observations wake formation are 
those given who employed stereoscopic 
cinematography record the fluid movement behind 
sphere being drawn steady motion through water with 
Reynolds Numbers 150 10,000. The vortex ring for- 
mation and shedding were obvious 450, and, these vor- 
tices moved askew giving the impression being linked 
together “vortex chain” (“Wirbelkette”). 
Reynolds Number about 1000 there occasionally ap- 
peared discrete vortices which with increasing values took 
the upper hand over the vortex chain until 1500. From this 
value upwards only discrete vortices—referred 
“vortex noted. These vortices—not 
confused with the Von Karman vortex street occurring 
the wake cylinder—were unstable, and they joined 
together form periodic balls vorticity, “Wirbel- 

The rhythmical and therefore non-turbulent compon- 
ent the cyclic vortex usually described 
the Strouhal Number Fd/U,, where the 
shedding, and Moller presented two Strouhal 
Number curves which are reproduced Figure Curve 
gives the variation frequency vortex shedding with 
increasing flow, while curve indicates the frequency 
formation the Winny reported 
Strouhal Numbers considerably lower than those given 
Moller, but his results not seem 
Schmiedel’s results for freely-falling spheres indicated 
Strouhal Numbers 0.327 and 0.389 690 and 955 
respectively, which are not excessively higher than 
curve, view the experimental difficulties. 
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Figure 3—Sphere wake configurations 
at, and above the upper critical 


Reynolds Number (Foch and Chartier). 


Foch and Chartier developed 
stereoscopic method for obtaining quantitative survey 
the fluid movement the wake sphere air 
current, the time-lapse photography tracer particles 
injected into flow. Figure and were obtained 
from the results for Reynolds Numbers below, at, and 
above the higher critical value. the subcritical region, 
the wake clearly consisted two intertwining vortex 
whose points separation from the sphere 
moved about the latter’s surface. The axis symmetry 
the wake spun about the sphere axis parallel the flow, 
the diagrams. Figure was obtained the higher critical 
value, and here the wake had inner “dead water” core 
demarcated the dotted line, which the tracer par- 
ticles moved coil-like direction transverse the wake 
axis symmetry. higher Reynolds Numbers, the flow 
probably representative sphere systems 
which there has been laminar-turbulent transition the 
attached boundary layer—seemed generate spiral wake 
given Figure 3c. 

surprising result was recently reported Cometta 
the course hot-wire investigation the 
unsteady flow pattern air the wake spheres 1.9, 
2.2, 3.7 and 5.7 cm. diameter. contradiction the 
smooth Sr-vs-Re curve obtained Cometta’s 
results showed abrupt change the Strouhal Number 
Reynolds Number 7,400. Although makes 
mention detailed results, some his limited 
data seem fall the curve given Moller for the 
“Wirbelkette” and others obtained with larger spheres 
near the curve for the “Wirbelvolken”. Cometta 
loss explain this apparent effect sphere size, but 
comparison with Moller’s curves suggests that may 
measuring the fluctuations different regions the 
wake the two cases. Cometta also the opinion that 
helical wake movement never appears, although this 
argument difficult uphold the light Foch and 
Chartier’s detailed measurements. 


Sphere Wake Separation Circle and Pressure 
Distribution 


The position the circle formed the separation 
the boundary layer from the sphere surface changes with 
the Rey Number and this phenomenon exerts 
dominant influence the vortex formation and the over- 
all momentum transfer problem. Taneda studied the 
450. Data obtained with different sphere diameters all fell 
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smooth logarithmic curve which showed the angle 
between the rear stagnation point and the separation 
72° 450. This subject was investigated the 
same time Garner whose technique implies that 
the boundary layer for mass transfer coincident with 
the velocity ‘boundary layer that the minimum rate 
solution soluble sphere marks out the separation 
circle. The results differ from those Taneda below 
100, Garner’s values being higher; but above 
100 the results are almost coincident. Garner’s data show 
remain constant above 500, but this contra- 
dicted earlier work Ermish who found, 
means pressure measurements the sphere surface 
the range 800 26,000, that increases steadily 
throughout the range. The discrepancy may due the 
deformation the dissolving sphere surface the vicinity 
the separation point. Ermish felt that the increase 
may account for the gradual and moderate increase 
this region. This drag coefficient increase may also 
caused shifting towards the sphere the laminar- 
turbulent transition point the separated boundary layer 
The normal pressure the sphere surface 
maximum the front stagnation point and falls smoothly 
minimum the separation point; then rises 
amount which inversely proportional the drag co- 
efficient, and then remains roughly constant across most 
the width the wake. The form drag may 
measured directly from the aphical integration the 
pressure distribution the vicinity the sphere 
and Flachsbart (23) used this method obtain value 
0.460 for the drag coefficient relatively high Reynolds 
Numbers, which agreed exactly with 
obtained from force measurements—thus illustrating the 
negligible role played frictional resistance 
region. 


Theoretical Studies the Sphere Wake 


The theoretical problems resulting from the separation, 
pressure distribution, and wake shedding phenomena are 
fundamental, and such complexity that they remain 
properly developed for systems far less complex than 
that the sphere. The application what has resulted 
from investigations simpler flow systems hampered 
insufficient experimental evidence the nature and 
characteristics the sphere ambient flow. 
made available present measuring techniques are 
cult interpret, being ambiguous and often contradictory. 
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Theoreticians have attempted provide crutch for the 
experimentalists who work the wake region pre- 
dicting from “reductio absurdum” reasoning what 
vortex existence would very difficult justify. 
The studies, some which were initiated Karman and 
Rubach deal with the different vortex 
head however, cautions that the mathematical form 
these disturbances may not correspond those which 
exist nature, and that the may thus quite 
infinitesimal and finite perturbations will discussed 
detail following paper dealing with the effects free 
stream Vorticity ). 

There are many possible types vortex systems 
which can investigated, but the accent has been 
which had seemed detected experi- 
ment. The most obvious system would one made 
series separated rings moving their axis 
symmetry. Levy and Forsdyke proved that such 
planes was only partially and inferred that their 
production occur only transition stage more 
stable systems. 


For Numbers the wake cylinder 
consists two stationary and symmetrical the 
range 10°, asymmetry sets and eddying 
alternatively from ‘side side, pproximated Karman’s 
“vortex street” which proved stable for 
one configuration only. The striking success achieved for 
the cylinder aroused hope that could extended 
somehow the sphere, which after all its axisymmetric 
analogue; Rosenhead has shown 
Karman vortex street three-dimensionally unstable when 
allowance made for the finite thickness the vortex 
filaments. Jeffreys points out that the 
slight when the vortices are thin, they are initially. The 
conclusion is, however, that there simple extension 
the Karman vortex street three dimensions. 


those 


Low 


Levy suggested quasi-steady motion, consisting 
two interlocking helical vortices, but Jeffreys 
showed that would very difficult maintain such 
The also erroneously con- 
cluded that the actual motion the wake was probably 
completely turbulent, i.e. there was periodicity asso- 
with it. Lan demonstrated the instability 


vortex system. 


pariello 


helical vortex filament, but Levy and Forsdyke 
countered showing that the helices are stable above 


Sranton 


against a 


Rosenhead pointed out theoretical reasons 
vortex proposed Levy, and indeed 
the the 
which had originally seemed indicate 


vortex showed closer inspection made 


helical 


photographs obtained former wake 


fiat pate 


from the 
there much occurring the sphere wake 


not presently discerned, alone understood. 


been 


that some the conflicting 
which the Reynolds Number 


phenomena occurring with different 


properucs when appl lied to fluids of greatly different 

that have been found sufficiently good for one 

inapplicable another” and interesting 


and gas are discussed recent provocative paper 
Evans 


There are indications that acoustic phenomena result- 
ing from the high-frequency oscillation shed vortices 
may produce beat frequencies which will dictate the value 
the Strouhal Number and through the value 
the drag coefficient. Experiment will have 
out the finer details the system which will include the 
fluctuating components the pressures and drag forces 
the sphere. 


The Two-Dimensional 
Observations 


The quantitative characteristics the flow perpen- 
dicular infinitely -long cylinder can times serve 
useful first approximation the flow about sphere, 
even though the field the latter not general axisym- 
metric analogue the former. The variation the flow 
patterns well the drag coefficients the cylinder 
with increasing Reynolds Numbers similar 
sphere and the lower and higher critical Reynolds Num- 
bers occur approximately and 105, respectively, 
Since the cylinder represents considerably simpler 
theoretical and experimental attention. Indeed some work 
done with cylinders may never repeated with spheres 
without the introduction inv alidating effects resulting 
from the perforation the flow field and the restraining 
action the supporting apparatus. The problem the 
cylinder therefore vital and immediate antecedent 
the more difficult sphere system. 

Studies the wake cylinder may grouped 
into two classes: the first deals with the wake the 
immediate vicinity the body, while the second con- 
cerned with the downstream development approxi- 
mately diameters more. The latter class, which has 
been the subject work Townsend Kobashi 
Grant and others, related, but not 
importance the problem the evaluation the drag 
coefficient, since the fluid pressure region not 
far removed from that the free stream. 

Papers given Ahlborn 1902 and Bernard 
1908 were among the first give photographic evidence 
periodic arrangement vortices the wake 
cylinder, although knowledge this had already been 
demonstrated Leonardo Vinci. Karman’s classical 
treatment the vortex street which was presented 
1911 initiated extensive series theoretical investi- 
gations the alternate system vortices which has yet 
failed relate properly ‘the downstream flow the wake 
the initial conditions the surface the body. The 
fascinating geometry the periodic flow and the theo- 
workers its study. The 
Reynolds Number relationship for cylinders —an effect 
which cannot accounted for the Karman theory 
since the latter assumes fluid—is much better 
established than for spheres, and this due part the 
work Carmichael, Dupin and Fage 
and Johansen Fage and Relf 
Relf and Owen Rosenhead and Schwabe 
and recently Etkin, 
and Keefe Roshko Hama 
Richardson and Cometta 


(50) 
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The flow the immediate vicinity the front 
the cylinder, ie. the front “stagnation point”, has 
been assumed free from fluctuation, and 
this seems confirmed all photographs the 
region, although some fluctuations can observed 
cinematographic studies. Piercy and Richardson stud- 
ied the flow very the surface cylinder 
means hot-wire anemometer, and their very precise 
underline the great complexity the flow. Surprisingly, 
there are large velocity fluctuations very close the 
front stagnation point, ‘and the disturbance seems 
damped out and then disappears into the boundary layer 
towards the wake. There has been very little 
comment the literature regarding this phenomenon, 
may important precursor the wake 
The results also indicate that the vicinity 
the separation point there arises “funnel vorticity” 
whose very large-amplitude fluctuations may have been 
fed the agglomeration smaller vibrations occurring 
the separating boundary layer. The funnel seems 
the starting point the large vortices the rear. 


The fate the boundary layer after separates from 
the cylinder has received better attention 
sphere analogue. Fage and Johansen recognized that 
the river vorticity emanating from the separation point 
and flowing into the heart the wake represented one 
the key phenomena which must incorporated into the 
understanding the overall momentum transfer. They 
investigated carefully the structure the separated layers 
springing from the surface cylinder and various other 
two-dimensional bodies measuring their velocity pro- 
files and then calculating their vorticity content from 
this information the followi ing manner: given point 
two-dimensional flow, vorticity given the 
expression: 


where and are the instantaneous velocities the 
mutually perpendicular and directions respectively. 
the flow the separated were free from fluctuations, 
then the total vorticity passing through its cross-section 
per unit time would given by: 


yi 
where and are respectively the inner and outer limits 
usually expressed dimensionless term obtained 
dividing the square the free stream velocity, 


K/U,? (3) 


the vorticity moves deeper into the wake acted 
viscous forces and turbulent and grad- 
ually annihilated vorticities opposite sign 
fore ata given point x downstream from the separation 
point the flow characterized the vorticity survival 
fraction defined 


(4) 


the subscripts and referring values and the 
point For the cylinder, Fage and 
Johansen found that 1.01 and that falls ste adily 
0.93 one-half the cylinder diameter distance 
from the separation point. has been suggested that 
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Figure 4—Velocity fluctuations vicinity cylinder sur- 
face and Richardson). 


the point formation the vortex centres, drops 
value approximately 0.5 but there are yet 
detailed measurements the variations this parameter 
with position the wake and which the flow Reynolds 
latter dictating the average ratio the 
turbulent vorticity the total. 

should noted that the foregoing treatment com- 
pletely ignores the vorticity associated with the non-ste ady 
component the flow. Fage and Johansen felt that this was 
justified that the velocity fluctuations the separated 
flow—which they me: asured with relatively large and 
uncompensated hot-wire and recorded with string gal- 
vanometer—were sufficiently allow (2) 
ization, since large scale fluctuations can pres- 
ent shown the work Pierey and Richardson 

results Schiller and Linke and 
show that the free boundary layer the cylinder will 
remain laminar distance which diminishes with 
increasing Reynolds Numbers and free stream turbulence, 
but the mechanism for this behavior obscure. 
understanding the separated boundary layer vital 
the solution the overall problem 
and long range program directed toward this end 
reported preliminary observations deal with 
the stability boundary separating from sudden 
that curvature, pressure distribution and free 
stream turbulence will investigated. Similar observations 
Maekawa and reiterate the need for the 
separated layer 
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The separated boundary layer and its dev elopment into 
the wake were studied Roshko detailed 
survey with hot-wire anemometer. Three regions were 
evident: Revnolds Number there was the start 
“stable” range which could discerned change 
the attached rear vortices 
forming the vortex street, but 
minar-turbulent instability present the 
enure rom 150 300 (called the transition 
range) there was Jaminar-turbulent instability, but this 
occurred the boundary the transition 
approaching the separation point with increasing Reynolds 
Above 300 the irregular range 
transition occurred the separation point, 
velocity fluctuations had both periodic and truly turbulent 
components, which would correspond the akes shown 
photographs for falling spheres Reynolds 
Roshko terms the fluctuation 
observed particular point being given 


Numbers above 1000 | 


from the truly turbulent fluctuations, 
the fluctuation contributed the shedding frequency 
the fluctuation resulting from the 


where results 


both and corresponds the frequency 
2N. The higher harmonics were found negligible. 
corresponding the above velocity The 
energy spectrum, the energy variation with frequency, 
shown the hot-wire anemometer exploration not 


smooth one the beginning the wake, and 


exhibited and There was, however, 
tendency for the spectrum smooth out 


transfer from the high 


noved Qownstrcall bY 


bands to the lower ones. 


reported asberg and Cooper 
have revealed spectra whic correspond shape 
those obtained The range 
characterized but non-turbulent com 
ponent the wake energy has been recently shown 
and corresponding Strouhal Numbers meas- 
will made the section dealing with 
the effect high turbulence inten 
gi ‘ t or whe Cparalé d boundary layer ol 
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stable stationary vortex region was qualitatively similar 
the results for the sphere given Figure they 
show that given Reynolds Number, the cylinder wake 
more than its sphere counterpart. 


Taneda found what termed remarkable fine scale 
phenomenon the region 45. Very small 
disturbances seemed appear the rear centre the 
moving downstream the interface the two 
stationary vortices, the end which they trembled for 
short time and died away. Reynolds Numbers above 
they passed downstream and may have seeded the 
periodic shedding. Although not mentioned Taneda, 
there some possibility that this phenomenon linked 
with the instability noticed just above the front stagnation 
point Piercy and Richardson 


Acoustic Phenomena 


There good reason believe that fine scale phe- 
nomena the flow which have yet fully dis- 
covered may play dominating determining 
the properties the wake. These are probably linked 
with the acoustic radiation which accompanies fluctuations 
velocity, the aerodynamic sound differentiated 
from that resulting from bodily vibrations. Richardson 
recently presented eloquent plea for reunion 
efforts dealing with turbulence and acoustic phenomena, 
emphasizing the aspects common both subjects; the 
importance the relationship the acoustic radiation 
resulting from isotropic turbulence the general problem 
interesting paper Meecham and Ford 

immediate importance the wake problem 
theory developed Lighthill 1952 which indi- 
cates that the velocity fluctuations turbulent wake 
and the fluctuating stresses the body surface act 
sources generating sound hose intensity varies 
exponentially with the relative fluid velocity. The ‘difficult 
measurement the acoustic intensities for different loca- 
tions the cylinder flow field has been the subject 
recent investigations Bingham, Weimer and Griffith (7), 
Gerrard Phillips Etkin, Korbacher and Keefe 
Richardson and McGregor has been found that 
fluctuating pressures caused the shedding vortices 
occur all along the cylinder surface, the intensities 
which rise from the front stagnation point maximum 
the sides. results, which include the evaluation 
the fluctuating component drag coefficient, are 
but yet incomplete. They inay also ex- 
pected explain the rather marked effects sound fields 
noted the mass and momentum transfer 
frequency blunt bodies. 


Theoretical Treatment the Cylinder Wake 


the evolution theoretical description 
the wake phenomena the cylinder yet 
state, the work that has been done valuable 
towards the eventual the sphere system. 
are three Conventional avenues approach the problem 
which are due Karman 


and Lees 


and Crocco 


notes that the flow about cylinder 
“ideal” fluid the potential flow, 
methods 


this would predict 


can readily calculated the 


however, 


the body the that the down- 
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Figure 5—A comparison experimental Strouhal Numbers with Shaw’s theory. 


stream side, that the net drag force which follows 
from the integration the distribution pressure the 
body surface would zero. Kirchhoff takes more 
realistic proach dividing the flow into two regions 
consisting outer potential flow field the wake 
envelope, and the wake region itself. The two regions are 
separated shear laver which the theory ‘idealized 
sheet discontinuity. The pressure distribu- 
tion calculated the wake boundary 
flow roughly approximates the observed 
one, but the theory runs into difficulty 
important wake region that assumes the velocity 
the fluid the free stream velocity, resulting 
base pressure the wake region which equal 
the free stream pressure. The predicted drag the body 
therefore much lower than the value. 

The celebrated Karman vortex street theory 
other hand attempts relate the base pressure with 
building and shedding vortices the wake. 
periodic wake Karman arrives the relationship 


Api 


where the drag the fluid density, 
spacing vortex rows igure the distance 
between consecutive vortices, UU. the mean stream ve 
locity. the the vortices relative 
velocity and the vorticity content 
ual vortex, (U, stability con 


siderations and assuming ideal fluid behavior 
deduced that can only have value 


asa result of hus momentui balance, it Was assuined that 


« Was piven by the CX pression tanh (gh a) so Chat iti 


could then found the 


Goldstein points out that even Karman’s two as- 


sumptions are valid, the drag could 
evaluated unless both and were experimentally 
determined. the event that the Strouhal Number were 
known, could evaluated from the expression 


and the vorticity per unit would equal 
K (U, U, 

and understandably what the author calls 
“disturbing” modifications to the von Karman theory have 
recently been proposed After citing evi 
dence that the Karman configuration h/a has 
higher ordes instability, shows that h/a ratios are 
stable providing that values and are used and 
that the instability present only with the 
periodic ity. bos NONVISCOUS sy steilis h/a Is proved ta be 

the same paper cites the possibility that the 
to the traisverse acceleranon of the wake fmuid, he 


the expression 


the Wake, and k thie Patio ol the louosel\ de fined ake 
diameter d* (usually taken as beimy equiv alent to h) to the 
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have assumed predict observed Strouhal Numbers 
are not unreasonable. 


Next follow attempts combine the Kirchhoff and 
Karman theories, not theoretically solve the drag 
problem least reduce the number parameters 
which must supplied experiment. Heisenberg 
assumed the Kirchhoff flow field what first appeared 
successful attempt predict the ratios and 
a/d needed for the Karman drag coefficient. The agree- 
ment his results with experiment have been shown 
fortuitous his derivation having greatly under- 
estimated the complexity the flow. 

points out that theory should have its wake 
fluid velocity corrected and Prandtl and Roshko 
have attempted this replacing with fictitious 


fluid velocity U,, defined by: 


C,, called the wake underpressure coefficient, related 
the drag coefficient and calculated from the differ- 
ence the mean fluid pressure the wake and the 
free stream pressure the relationship 


used Equation (11) very approximate 
theoretical derivation showing that 


and Birkhoff incorporated the vorticity survival fraction 
into this equation that has the final form: 


This equation would undoubtedly have refined 
before could considered theoretically signifi- 
cant. 
Birkhoff added the above concepts the Karman 
treatment, noting that the Karman formula for can 
rewritten give 


where the ratio needed for the Karman 
formula. substituting for Equation (14) this leads 

This expression cannot yet confirmed experimentally 
since adequate measurements are lacking. Birkhoff 
also introduced modified Strouhal Number which 
uses the actual wake diameter place the body 
diameter 


Sr, Nd*/U, (16) 

assuming and using the relationship 
Sr, = 6h (1 €)/a (17, 


(14) ending with the expression 


(Bk?) 


proposes modified Strouhal Number which 
incorporates not only the wake diameter but also the 
ficticious wake velocity 


and derives expression for which effect 
identical derivation and result Birkhoff’s Equation 
(17) except that replaced alone rather than 
Bh. then employs Equations (17), (18) and (19) 
graphical solution which amounts compact summary 
the results combining the Kirchhoff and Karman 
theories. substituting the given (18) into the 
Karman drag Equation (7) family curves having 
d/h the ordinate, the abscissa and the vary- 
ing parameter are obtained. Next, using the empirically 
modified Kirchhoff theory i.e. calculating the flow 
the wake the methods classical hydrodynamics and 
arbitrarily substituting different values the base 
pressure reflected and then calculating the drag 
coefficient from integration the pressure about the 
body, Roshko obtains family curves having d/d* 
the ordinate, the abscissa and the shape the 
cylinder the varying parameter, thus extending the 
treatment cylinders non-circular cross-section. 
the assumption allowed that then the points 
intersection the d/h vs. and d/d* vs. curves 
represent wake solutions for varying values and what 
most interesting, varying cylinder shape. This means 
that only need measured determine d/d* and 
for any cylinder shape, allowing the drag coefficient 
determined from plot vs. obtained from 
the modified Kirchhoff theory. Roshko points out that 
any interference the wake the form other bodies 
free stream turbulence would change and therefore 
and illustrates this with some initial experimental 
evidence. 

The times superficial treatment given the wake 
Birkhoff and Roshko seems have exhausted that 
which can done with the limited experimental results 
presently available. Further empirical work now pro- 
gress should allow replacement some the parameters 
with others which correspond more truthfully 
the actual wake phenomena. 

third theory which may useful, but 
not yet been dev eloped for subsonic wakes blunt 
bodies that momentum diffusion involved 
mixing theory for the interaction between dissipative flows 
and nearly isoentropic streams proposed Crocco and 
Lees considers the difference pressure between 
the base the wake and the free stream resulting from 
the diffusion momentum across the discontinuity 
velocity forming the wake envelope. cannot 
successful since omits vital part the mechanism, 
demonstrated Roshko, namely that the forma- 
tion the vortices, but can probably incorporated 
into overall solution the wake. 


The theories thus far presented are based conven- 
tional fluid flow concepts, and entirely possible that 
the difficulties encountered are part due omission 
less evident but vitally important occurrences the 
wake. theory such nature was presented Shaw 
been avoided overlooked nearly all the many 
workers the field who have published since—this spite 


its apparent success providing completely theo- 
retical derivation the sought after ratio 
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essence, Shaw argues that there are three disturbance 
centres the flow about cylinder, two strong ones 
the separation points and third weak one the front 
stagnation point—a model which complete agreement 
with Piercy and Richardson’s results (Figure These 
centres emit acoustic pressure pulses which travel the 
cylinder boundary layer and through the ambient flow 
sonic velocities. From the time taken the pulses 
travel between disturbance points, Shaw shows their fre- 
quency about equal where the speed 
sound. Moreover since shed vortex continues 
oscillate this frequency sends back pulses the still 
attached layer frequency which due the Doppler 
effect reduced the ratio (U, U,)/C] thereby 
interacting with the acoustic vibrations the attached fluid 
vibrating with the original frequency, and consequence 
producing strong pulses with the beat frequency 
which turn cause the shedding the 
vortex. The shedding frequency vortex pairs there- 
fore that the Strouhal Num- 
ber given 


This effect saying that 


result that agrees very well with experiment and also 
with Birkhoff’s proof that “a” given wake invariant. 
Figure shows the comparison the theoretically pre- 
dicted Strouhal having been obtained 
inserting experimental values into the Karman equation 
some the available experimental results and the 
theoretical curve almost appears have been drawn 
the method least squares. Excellent results have also 
been obtained for cylinders non-circular cross section 
revising the passage time for the pulses result 
the shape the body. Comment and criticism 
Shaw’s exciting theory will awaited, particularly from 
those active the the aerodynamic sound 
generated the environment immersed body. 


summary, review the sphere wake literature 
shows that the comprehension this complex region 
impeded lack adequate experimental substance 
combination with the theoretical difficulties common 
studies flow systems which the inertia terms are 
significant. There immediate need for paralleling 
application the sphere wake the albeit limited and 
unrefined information obtained from the simpler two- 
dimensional cylinder system. spite its shortcomings, 
the information presented forms foundation for the 
understanding the aerodynamics particles moving 
fluid environments where complications 
which include the particle acceleration, shape, rotation, 
concentration, and the fluid turbulence become operative, 
and these matters will dealt with the forthcoming 
installments this series. 


Nomenclature 


Any consistent set units may employed. listed 
are merely illustrative. 


distance between consecutive vortices, 

body diameter, ft. 

wake energy associated with fundamental shedding 
frequency, 

wake energy associated with first 


shedding frequency, ft. 


of 25) Rosenhead, 


wake energy associated with truly turbulent fluc- 
tuations, ft. 

shedding frequency sphere wake, 

lateral spacing between vortex centres for cylinder 
wake, ft. 

vorticity passing through separated layer cross- 
section per unit time, Equation (2) 

ratio ficticious wake velocity fluid velocity, 
dimensionless. 

ratio wake diameter cylinder diameter, dimen- 
sionless. 

average length wavering portion wake, 

from rear stagnation point vortex 
centre, ft. 

length sphere wake, ft. 

shedding frequency cvlinder wake, 

first harmonic shedding frequency, 

mean static pressure free stream, 

mean static pressure wake flow, 

fluid drag force body, force, (for 
two-dimentional case 

Reynolds parameter, dimensionless. 

Strouhal Number, dimensionless, Figure 

Sry Strouhal Number due Birkhoff, Equation (16). 

Strouhal Number due Roshko, Equation (19). 

free stream velocity, ft./sec. 

velocity vortices relative stream velocity, 
ft. 

ficticious wake velocity, ft./sec., Equation (10). 

total fluctuating component fluid velocity 
direction, ft./sec. 


tion associated with the fundamental 
harmonic the shedding frequency, and the truly 

co-ordinate parallel stream direction. 

co-ordinate transverse stream direction. 

angle between rear stagnation point and separation 
circle, degrees. 

vorticity survival fraction, dimensionless, Equation 

(4). 

vorticity survival fraction point dimensionless. 

half-angle wake cone, degrees, Figure 

ratio dimensionless. 

vorticity content per vortex, 

reduced vorticity, dimensionless, Equation (3). 

reduced vorticity separation point, dimensionless 

reduced vorticity point dimensionless. 
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Development Design Criteria for 


Biological Treatment Industrial 
Effluent Containing Waste 


MILLS? 


The development design criteria for waste 
water treatment plant described. Operating results 
for laboratory activated sludge unit, Dowpac 
trickling filter pilot plant, and aero-accelator acti- 
vated sludge pilot plant are given. Studies with the 
aero-accelator pilot unit indicated that observing 
certain critical nutrient and dilution requirements 
the industrial effluent, containing toxic 2,4-D waste 
water, could satisfactorily treated bio-oxidation. 
phenol removal 99% and biochemical oxygen 
demand removal 90% 95% was ob- 
tained. the full scale treatment plant, the esti- 
mated cost operation was calculated $0.50 per 
pound phenol removed, and $0.08 per pound 
removed. 


AUGATUCK Chemicals, Division Dominion Rubber 

Co. Limited, manufactures wide variety rubber, 
agricultural, industrial chemicals its plant 
Elmira, Ontario. The industrial waste 
operations currently averages about 0.238 million imp. 
gallons per day and has five- day biochemical oxygen 
demand averaging about 700 p.p.m. 

The “phenolic” wastes are presently chlorinated and 
discharged effluent ponds having overflows, while 
the non-phenolic wastes are discharged equalizing 
lagoons. 

This discussion intended outline investigation 
method improve present treatment facilities the 
industrial wastes shown Table The investigation 
started the laboratory and advanced the plant 
phase. Two different biological pilot units were used 
during the study; the trickling filter and activated sludge 
system. The laboratory and pilot plant trickling filter 
studies include the treatment the phenolic waste waters 
only. The pilot plant activated sludge studies include the 
treatment the principal component waste streams shown 


Table 


The history the approach this problem dates from 
about 1950. Among the various types industrial organic 


1Manuscript received July 6, 1959. 

2Naugatuck Chemicals Division of Dominion Rubber Co. Limited. 

paper presented the 14th Purdue Industrial Waste Con- 
ference, Poke University, Lafayette, Indiana, U.S.A. May, 1959. 
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wastes produced, the most troublesome and costly from 
the standpoint treatment, dichlorophenol 
containing waste resulting from the production 
dichlorophenoxyacetic acid (2,4-D). Several different 
methods treating this waste water have been investi- 
gated attempt improve treatment efficiency and 
cost operation. 


The alkaline chlorination process was adopted about 
1951 and still operation today achieving 95-98% 
destruction the dichlorophenol content the waste. 
Although the degree treatment considered satisfac- 
tory for the chlorination process, the treated waste still 
contains 20-30 p.p.m. dichlorophenol and hence cannot 
discharged the receiving stream. The cost treat- 
ment about $1.00 per pound dichlorophenol destroyed. 


Laboratory Investigation 

1954, was decided investigate the bio-oxidation 
process cheaper and more efficient method treating 
the “phenolic” waste waters. “Phenolic” wastes are in- 
tended here include both the phenol and 2,4-D com- 
ponent waste streams shown Table 
studies began with the investigation synthetic waste 
waters containing phenol, and/or 2,4-D 
acid. this, the study advanced investi- 
gation the actual plant phenolic waste streams. The 
laboratory units, Figure were designed permit con- 
tinuous feeding synthetic and plant wastes with con- 
tinuous removal the treated effluent. Each unit was 
seeded with sewage activated sludge and operated with 
continuous feeding and aeration, hours day, days 
week. Feeds were changed, samples for analy sis taken, 
temperatures recorded, and any changes operating pro- 
cedure made the same time each day. Feed waters and 
were chemical oxygen demand 
and where possible the individual contaminating 
component the waste. Nitrogen ammonium nitrate 
and phosphorous dibasic phosphate were 
added the feed give approximate C:N ratio 20:1 
and C:P ratio 

Based the laboratory investigation both 
and actual plant ‘phenolic” wastes, was found that, like 
phenol, dichlorophenol and 2,4-D acid could 
one the laboratory units treating waste water 
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949), 
Res, 


TABLE 


Phenol 
gal 
min. p.p.m. p.p.m. Ib.’ 
400* 63* 1330 210 
Misc. chemicals 100 350 500 
Polyester resins 1020 200 
*Phenol (Carbolic id). 
(a) Investigation 2,4-D Waste 
The pilot unit was seeded with sewage activated sludge. 
shown Table Based the encouraging laboratory Using the phenol waste water influent, the biologically 
results, was decided advance the pilot plant stage active slime was established the packing and the 
the investigation. activated sludge tank. After activation, the influent was 


gradually changed the dichlorophenol containing 2,4-D 


Pilot Plant Trickling Filter Investigation waste water, which had been added nitrogen and 


The pilot plant trickling Figures and was phosphorous salts supply minimum carbon-to-nitro- 
designed combined trickling filter and activated sludge gen-to-phosphorous ratio During month 
system employing Dowpac 10, formerly Dowpac period attempt was made treat 2,4-D waste water 
the filter media the unit, but without the desired success. 
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Figure diagram for pilot plant trickling filter. 
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TABLE 


OPERATING SUMMARY LABORATORY BIOLOGICAL UNIT 
2,4-D WASTE WATER 


Influent Effluent 


Flow, influent, litres/day 


Flow, dilution, litres/day 


oo 
w 
Co 
w 


D.C.P., p.p.m. (after dilution) 148 138 


co 


p.p.m. (after 1045 970 240 195 


removed cu. ft./day 175 150 


When operating 2,4-D waste alone, suspended solids synthetic waste water containing only dichlorophenol. 
could not held the activated sludge tank. The bio- The result was the same when treating the actual dich- 
logical growth the packing changed from gelatinous lorophenol containing waste water. This study 
material characteristic phenol waste treatment showed that dichlorophenol, although can classified 
rather hard pebbly slime. When any the biological compound can not satisfactorily oxidized 
growth the packing sloughed off, was not replaced. biologically under the same conditions 
one period operation the influent consisted The reaction the bacteria the 


cated also that the latter protoplasmic poison that 
inhibits some way the normal metabolic processes 
the bacteria. Chlorine itself well known bactericide 
and there was evidence here indicate that the introduc- 
tion chlorine into the phenol molecule creates com- 
pound, dichlorophenol, that apparently toxic the 
bacterial culture which existed during the study. 


Actually the pilot unit treated 400-500 gallons 2,4-D 
waste water per day. Before being applied the filter, 
however the waste was diluted with water. During the 
most efficient period operation, average 7.8 per 
86%. addition, average 16.7 Ib. per day 
chemical oxygen demand were removed show- 
ing average removal 70%. summary the oper- 
ation the unit treating waste water shown 
Although operation the pilot unit 
waste could not called failure, could not con- 
sidered commercially feasible, since the unit could neither 
treat sufficiently large volume the waste, nor could 
the dichlorophenol content the reduced 
below that which could achieved 

the end the production season, portion 
the phenol (carbolic acid) waste water shown 
was very successfully treated the unit. average 
11.8 per day phenol were removed showing aver- 
age removal addition, erage 30.3 Ib. 
per day was removed showing average re- 
moval summary this operation shown 


Figure 3—Pilot plant trickling filter. Although was not the ultimate objective, 
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TABLE 


OPERATING SUMMARY PLANT (Dowpac 10) TRICKLING FILTER 
2,4-D WATER 


Influent Packing Effluent Aeration Tank Effluent 
Total 
Percent 
Max. Av. Max. Av. Av. Max. Av. Av. 
Removal Removal 
Flow, dilution, gal./min. 4.0 3.0 4.0 3.0 4.0 3.0 
p.p.m. (after dilution) 780 585 215 170 215 170 
Max. 
TABLE 
OPERATING SUMMARY PLANT (Dowpac 10) TRICKLING FILTER 
PHENOL WASTE WATER 
Influent Packing Effluent Aeration Tank Effluent 
Total 
Percent 
Max. Av. Max. Av. Max. Av. 
Flow, dilution, gal./min. 1.4 1.0 1.4 1.0 1.4 1.0 
Max. Av. 


| 
| 
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Figure 4—Aero-accelator pilot plant—17 cu. ft. aeration 
volume. 


was estimated that increasing the packing capacity 
factor adding three additional towers, the 
total volume the phenol waste water could treated 
the Dowpac pilot unit. The estimated cost treating 
the total phenol waste stream expanded unit was 
calculated $0.26 (summer operation) $0.34 (winter 
operation) per pound phenol removed. This cost in- 
cludes amortization, labor, power, steam, laboratory ser- 
vices and nutrients. 


Pilot Plant Activated Sludge Investigation 

Following the pilot plant trickling filter studies, 
became apparent that improved treatment facilities were 
required, not only for the phenolic wastes, but also for 
the non-phenolic wastes presently discharged the equal- 
izing lagoons. was therefore decided further the 
study investigating the activated sludge process. 

the general field industrial waste treatment there 
increasing recognition the advantages activated 
sludge treatment completely mixed system. Such 
system can, course, handle shock loads effectively, be- 
cause feed the system dispersed rapidly throughout 
the entire mass activated sludge. Furthermore, 
possible operate completely mixed system much 
higher loads than have characterized the traditional acti- 
vated sludge flow sheet the past. 


November 1957 further studies were undertaken 
with Aero-Accelator pilot plant. The pilot unit, 
Figure provided completely mixed aeration volume 
and clarification zone single basin. was in. 
diameter and had aeration volume approximately 
cu. ft. 


complete study program was set and carried out 
under the direction the pilot plant manufacturer 
The pilot work was complex problem involving the 
control temperature, air requirements, oxidation reduc- 
tion potential, dissolved oxygen, nutrient addition, degree 
neutralization, aeration solids, biochemical oxygen de- 
mand (5-day, 68°F.), and numerous other less important 
factors. The chemical and biochemical analyses carried 
out during this investigation were essentially those out- 
lined the 10th Edition “Standard Methods” 


Since previous studies indicated that the dichloro- 
phenol containing 2,4-D waste water would not itself 
support biological floc, was decided attempt treat 
this waste the Aero-Accelator pilot unit along with 
other more amenable waste waters. 


The laboratory and pilot plant trickling filter studies 
previously discussed, deal only with the treatment the 
phenolic waste waters. The pilot plant activated sludge 
studies discussed the remainder this paper, include 
investigation the treatment the principal com- 
ponent waste streams shown Table For this reason 
there direct comparison between the laboratory 
and trickling filter, and the activated sludge studies. 


Since the first four component waste streams shown 
Table would make about 80% the total organic 
and hydraulic loading full scale plant, pilot studies 
were limited these four waste streams only. The waste 
stream from the production polyester resins was not 
included the study since this waste did not exist the 
time the pilot plant investigation. addition, the 
time-honored method investigating 
waste stream was largely abandoned attempt 
shorten the study period. 


The Aero-Accelator pilot plant was seeded with sew- 
age activated sludge and satisfactory biological floc 
established influent consisting the aniline and 
nitrobenzol, and phenol waste waters. After activation, 
the proportion 2,4-D waste water the feed mixture 
was gradually increased until the aniline and nitrobenzol, 
phenol and 2,4-D waste streams, which nitrogen and 
phosphorous additions were made, were fed the unit 
the ratio their respective flow rates. The Aero- 
Accelator pilot unit operated for about two weeks this 
influent, and then for apparent reason, the character- 
istics the mixed liquor solids the aeration zone 
deteriorated and the floc was lost. The unit was again 
seeded and operated under the same conditions before, 
but again the floc was lost after two week period 
operation. Trace amounts cobalt, manganese and fer- 
rous salts were added the unit and this addition ap- 
peared have beneficial effect the characteristics 
the mixed liquor solids. this addition the initial period 
operation was extended and the Aero-Accelator pilot 
unit operated for another five weeks before any further 
sign floc deterioration was noticed. When this occur- 
red the trace element addition was replaced the daily 
addition settled sewage. The affect the sewage 
the sludge characteristics shown Figure where the 
settling rate plotted both before and after the sewage 
addition. The mixed liquor suspended solids concentration 
was essentially the same for both curves. The sludge 
settling rates were observed 1000 ml. glass graduate. 
The pilot unit continued operate for additional three 
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weeks the influent composed the aniline and nitro- 
benzol, phenol, and 2,4-D waste waters. Operation was 
successful during this period due the daily addition 
settled sewage 

Once the problem essential nutrient requirement 
was resolved, the study was able proceed and during 
the next phase the investigation the miscellaneous 
chemical waste stream was included the Aero-Accelator 
pilot plant influent. This waste was previously excluded, 
since was unknown from biological standpoint 
because previous laboratory pilot plant work had 
been carried out it. The influent first investigated 
consisted the miscellaneous chemicals, phenol, and 
aniline and nitrobenzol waste streams. Operation was 
excellent. The second influent investigated consisted 
the above three waste streams, plus the 2,4-D waste. Oper- 
ation was still good, but conditions operation were 
slightly more critical. summary the operation the 
Aero-Accelator pilot unit these two influents shown 
Table 

From the data obtained this last phase the study, 
was felt that the major characteristics the activated 
sludge system were known and full plant scale design 
could considered. 

The seven main factors which govern the design 
full scale treatment plant were calculated and are 
follows: 

The optimum temperature the operation was found 
range between 73° and 88°F. The median temper- 
ature was 85°F. 


The B.O.D. loading the activated sludge system was 
found vary considerably with the various raw waste 
streams. the typical waste-miscellaneous chemicals, 
phenol, aniline and nitrobenzol, and 2,4-D, maximum 
B.O.D. loading 375 per 1,000 cu. ft. aeration 
volume per day was obtained. Without the 2,4-D 
waste, maximum loading 790 Ib. B.O.D. per 1,000 
cu. ft. was possible, but this loading produced excessive 
sludge bulking. successful loading (no bulking) was 
200 Ib. B.O.D. per 1,000 cu. ft. While this loading 
somewhat higher than the average value shown under 


that lies within the range satisfactory operation. 


The air requirements for oxidation the organics 
with the bacteria operating the logarithmic growth 
phase, again when treating the typical waste including 
2,4-D, was 1.3 lb. oxygen per Ib. B.O.D. This oxy- 
gen requirement based substrate temperature 
and minimum mixed liquor dissolved 
oxygen content 2.0 p.p.m. 

The hydraulic loading the clarification unit pro- 
duce acceptable quality effluent from the Aero- 
Accelator pilot plant was found 0.24 gal./(min.) 
(sq. ft.) The controlling loading the Aero-Accel- 
ator pilot unit was, however, the B.O.D. loading, rather 
than clarification. The design overflow rate was select- 
from the sludge subsidence tests, rather than the 
hydraulic throughput for the pilot unit. believed 
that slightly higher rise rate can tolerated the 
full scale unit. this connection, interest 
note that the 2,4-D waste results considerably 
more dense and more rapidly settling sludge than 
produced when this waste not present the stream 
being investigated. 

Based total waste flow 165 gal./min., dilution 
factor 3.4 was found necessary for the process 
proceed successfully. 
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Figure 5—Sludge settling curve. 


The levels essential nutrients and essential metabolite 
supplement necessary insure active growth the 
biological system are follows: 

(a) The B.O.D.:N ratio must greater than 20:1. 
Lesser ratios may tolerated. 

(b) The B.O.D.:P ratio must greater than 100:1, 
and may vary 

(c) Essential metabolites the form sewage should 


estimated that the phenol content the treated 
will less than p.p.m. when influent con- 
centration ranges 500 p.p.m., less than 1.5 p.p.m. 
when influent ranges 1,000 p.p.m., and less than 
p.p.m. when the influent contains 3,000 p.p.m. 
phenol. B.O.D. reduction estimated the 
order 90-95%. 

Figure shows the principle dimensions the Aero- 

Accelator Type unit which was recommended for this 


Erelvent 
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Figure 6—Proposed aero-accelator activated sludge plant. 
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TABLE 


OPERATING SUMMARY PLANT AERO ACCELATOR 
ASTE WATER 


ION. = —~ — 
Influent No. Influent No. 
Effl. Effi. 
Removal 
Max. Av. Max. Av. Max. Max. Av. 
(50% effluent) 
B.O.D. (5-day, 68°F.) p.p.m. (after dilution) 860 730 375 330 
B.O.D., removed, cu. ft./day 790 350 375 130 
Influent No. Influent No. 
Max. Av. Max. Av. 
Clarification area overflow rate, 0.10 0.06 0.11 0.06 
Mixed liquor suspended solids, 8492 7566 8830 8075 
Influent No. Misc. chem. 100; Phenol 11; Aniline and nitrobenzol 
Influent No. Misc. chem. 100; Phenol 11; Aniline and nitrobenzol 2,4-D 
application. Some extra capacity was included provide 99.5% and chemical oxygen demand (C.O.D.) 
for overloads which may inadvertently occur. Basis for removal 95.3% was obtained. 
the design this unit, course, stems from the analysis Dichlorophenol, although phenolic compound, cannot 
of the pilot plant study. The operational characteristics be satisfactorily oxidized by bacterial action under the 
the unit are follows: same conditions can phenol. However, treating 
Flow capacity gal./min. the dichlorophenol containing 2,4-D waste water along 
ENT 
Hydraulic loading 0.3 /(min.) (sq. ft.) with other more amenable waste streams, the phenolic 
capacity content the waste can satisfactorily oxidized, 
loading—202 lb. cu. ft. aeration volume/day providing certain critical nutrient and dilution require- 
Nominal air required 1300 stand. cu. ments are observed. 
Maximum air 1350 stand. cu. preliminary assessment the operation the full 
Power drive scale Aero-Accelator waste treatment plant, was cal- 
The 1300 stand. cu. ft. air supplying culated that the total industrial including the 
3580 oxygen per day, and based 12% oxygen waste water, could treated for the same total 
transfer efficiency. yearly cost currently experienced for the treat- 
ment the 2,4-D waste alone. 
Conclusion 
From the laboratory and pilot plant investigations References 
described, design criteria was ultimately dev eloped for (1) Bryan, E. H., Molded Polystyrene Media for Trickling Filters, Indus- 
trial Wastes, 1, 2, 80 (Nov.-Dec. 1955). 
industrial effluent containing (2) Plastics Technical Service Bulletin, “Dowpac NF-90 and Dowpac 
toxic 2,4-D waste water. The ing significant facts (3) Technical (1954). 
, ; is i (4) “Stands Methods f the Examinati Water, -wage, an 
were established during this stu y Industrial Wastes”, 10 Edition, "APHA; AWWA; FSIWA (1955). 


The trickling filter design, employing Dowpac 
the filter media, satisfactory device for oxidizing 
phenol (carbolic acid). average phenol removal 
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5) Oginsky, E. L., Umbrait, W. W., “An Introduction to Bacterial 


Physiology”, W. H. Freeman & Company, San Francisco, (1954) 
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Detailed studies the air and flow patterns 
spouted wheat bed have been made column 
sizes ft. diameter. Orifice size, cone angle, 
bed height and air flow were varied establish 
their effect these flow patterns. Air flow within 
the bed was obtained from pressure drop gradients, 
and solid flow was obtained from particle velocity 
measurements taken the wall. Pressure drop, air 
flow distribution between spout and annulus, solid 
circulation rate and flow line data are summarized. 
correlation proposed which relates pressure drop 
and particle velocity the top the bed with super- 
ficial air velocity air flows above minimum spout- 
ing. Particle velocity measurements the spout were 
made from photographs, and were calculated approxi- 
mately from the drag equation for single particle. 


contacting fluids and granular solids spouted 
bed has received special study these laboratories 
recent years. Early work Mathur and de- 
scribed the mechanism spouting, and produced 
empirical correlation for predicting the minimum air 
required for spouting variety materials. Publications 
have appeared heat transfer (2) and potential applica- 
tions the drying wheat and wood chips Recent 
studies have been made the flow patterns air and solid 
within the spouted bed with the object developing 
rational correlations for predicting spouted bed behavior. 


The general flow wheat particles spouted bed 
illustrated Figure Wheat particles are lifted 
dilute suspension through the central core, high 
velocity air stream. the top the bed, where the 
central stream velocity suddenly reduced, the wheat 
particles fall back onto the bed. These then flow, 
slowly moving packed bed, downward and toward the 
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spout. The air enters the column through orifice, 
flows the spout and through the void space the 
moving packed bed. 

This paper summarizes the observations made 
column sizes ft. Pressure drop, minimum air 
flow, and the effect design variables air and solid 
flow patterns were measured, the design variables studied 
being bed depth, inlet diameter, cone angle and column size. 


EXPERIMENTAL 


Data were obtained from three columns: 0.5 ft. 
ft. glass, ft. ft. steel column fitted with plexi- 
glass windows and ft. ft. plexiglass column: 
0.5 ft. and ft. semicircular columns with plexiglass 
fronts were also used. Cone angles and orifice sizes were 
changed insertions between flanges provided for this 
purpose. Straightening vanes were placed the approach 
pipe the The air was filtered and then metered 
with calibrated orifice. 

Air and solid flow rates within the spouted bed were 
measured except that the static pressure 
measurements were made with Betz Manometer sensitive 
0.01 mm. water. 


Pressure Drop Across Spouted Bed 

The relationship between pressure drop and air flow 
shown Figure for wheat spouting ft. LD. 
column, in, and 1/2 in. bed depths. The section, 
the curve corresponds the pressure drop through 
fixed bed additional pressure required form 
internal spout within the bed. Wheat particles must 
pushed upward against the weight the bed. Void 
volume the packed bed surrounding the cavity re- 
duced and pressure drop increased. the internal 
spout moves the cone, point reached where the 
pressure drop decreases, (section BC), due the decreas- 
ing length packed bed between the internal spout and 
the top the bed. Further increases flow result 
sharp drop pressure drop where the internal spout 
breaks through the surface the bed. Point defined 
incipient spouting. flows beyond pressure drop 
remains approximately constant. 

When the air flow reduced, pressure drop follows 
the path EE’, E’, the spout suddenly collapses. 
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the cone just above the orifi pressure was noted 
cone height), the pressure drop ined 
the bed. cone and the top 
Attempt 
g 
pressure drop across the 
yout into the bed 
the pressure drop data Table low bed depths 
with average value 0.4 
gure 1—Spouted wheat bed. bed drop across unrestrained 
equal the weight the bed 
The air flov 
reproduced reliably. The section across spouted bed therefore 
drop flow relationship for the the pressure due the weight the 
equation for packed beds the pressure drop across bed 
hysteresis loop, and this occurs just abov 
the difference betw and show that spoutabl just above the maxi 
ased. The maxi- (1) the Ergun equation 
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TABLE 


SUMMARY PRESSURE AND SOLID FLow DaTA 
ARYING ATES AND BED IGHTS 


ft. 
3.13 47.3 1.75 0.074 165 
4.23 66.2 2.66 0.100 283 
| | | 
| | | | 
giving 


The maximum spoutable bed depth can now calculated 
equating the value for U,,, obtained solving the 
above quadratic, the standard spouted bed Equation 
(3). Figure shows the agreement between calculated 
and experimental values for variety materials reported 
earlier 
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1.0 3.0 4.0 


CALCULATED- FT. 


Figure 3—Measured vs. maximum spoutable bed 
depth. 


Minimum Spouting Velocity 


resulted empirical equation for minimum spouting 
velocity 


1/2 


Minimum air flows were measured for the ft. column 
using three different orifices and three different cone 
angles. The data are summarized Table For varying 
orifices the air flow required for minimum spouting in- 
creased slightly with orifice diameter. two fold change 
D,, (1-7/8 in. in.), caused about 10% increase 
minimum air required for spouting. This change 
agreement with previous data 

The effect cone angle minimum spouting velo- 
city not large. The minimum spouting velocity was 
similar for 45° and 60° cones but about 10% higher with 
the 85° cone. This difference larger than the experi- 
mental error, but not sufficient justify any change 
Equation (3). The calculated values based Equation 
(3) are shown Table These are about 10% lower 
than the observed values. 


TABLE 


INCIPIENT AIR VELOCITIES FOR 


3.0 2.60 1.55 
5.0 2.05 1.99 
4.0 1.88 1.78 
6.1 2.48 2.20 
2.33 2.03 
4.1 2.17 1.80 
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Figure 4—Effect column diameter air flow pattern. 


AIR FLOW PATTERN 


Air Flow the Annulus 


The distribution air between spout and annulus 
important assessing the effectiveness the gas-solid 
contact. Previous work with the 0.5 ft. column 
has shown that 70% the total gas flowing passes 
through the annulus. The distribution ft. column 
was therefore examined under various conditions. 

Air flow was estimated from pressure drop measure- 
ments different bed levels. relationship between 
pressure gradient and air flow was first established empiric- 
ally using packed bed. Measurements were also made 
pressure drop per foot bed height moving 
packed bed which air and solid flow was maintained 
the same value the annulus spouted bed. The 
pressure drop per foot vs. air flow curves were similar 
magnitude and shape both cases. Pressure drop profiles 
measured the annulus spouted bed were converted 
velocities using these curves and hence volume flows 
were obtained. 

Radial pressure gradients were small compared with 
the vertical pressure gradient that the horizontal com- 
ponent flow was neglected. The error introduced 
neglecting horizontal air flow significant only when 
L/r, ratios are small (L/r, 1). also assumed that the 
void fraction the moving annular bed the same the 
void fraction the stationary packed bed. This rea- 
sonable assumption the main body the bed, but 
not justified the cone region, the top the bed. 
the cone particles tend preferred orientation 
they are deflected into the spout. the top the bed, 
near the spout, the bed expanded and higher voidage 
results. 


Effect Design Variables Air Flow Pattern 


Air distribution within the annulus 
under varying conditions different columns. 
Figure shows the percentage the total air flowing 
through the annulus various fractional bed levels. The 
annular air flow increases with column diameter. Flow 
data are not recorded for the cone region the top 
the bed because the changes void volume these 
regions. 

Table the ratio Q,/Q, for the ft. column, 
tabulated several fractional bed the first 
group data, the cone angle and orifice size are constant, 
the bed depth varying from 3-6 ft. bed height 
increased, air flow required for spouting also increases, 
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but the same level bed approximately the same 
volume air flows through the annulus section. the 
3.3 ft. level the volume air flowing through the annulus 
and ft. beds was 257, 226, and 218 c.f.m. respec- 
tively. the second group data, bed height, and cone 
angle are constant, and orifice size varies. The ratio 
Q,/Q increases with increasing orifice size, particularly 
the lower half the bed. the third group, bed 
height and orifice size are constant, but cone angle varied. 
These data show air flows into the annulus 
particularly the lower section, when the cone angle 
small. the last group data, air flow was increased 
from 110 120% and 130% minimum spouting flow. 
increased, the annular flow remained essen- 
tially constant. the 4.2 ft. level for example, the 
volumetric flows were 252, 252 and 245 c.f.m. respectively, 
indicating that the excess air passed through the spout. 
This data annular air flow suggests that the most 
effective contact between gas and solids occurs when 
column diameter and orifice diameter are maximum, 
and total air flow and cone angle are minimum. 


Air Flow Through Spout 


Air flow through the spout different bed levels 
obtained from the difference between the total air flow 
and that calculated flowing through the annulus. The 
linear velocity obtained from knowledge the spout 
diameter. The effect design variables the air flow 
the spout may therefore ‘deduced from the preceding 
discussion annular air flow. 


Spout Diameter 


The spout diameter was measured from photographs 
the vertical section through the spout semi circular 
column. assumed that the spout was not significantly 
affected the wall and thus that any changes spout 
shape size were directly related changes design 

variables. Data are summarized Table where spout 
diameters are tabulated several bed levels for variety 
conditions. These data were taken from smooth curves 
drawn through the experimental points. 

Several general observations may made. For the 
column with in., the spout width decreased 
about 2.5 in. immediately above the orifice, then flared 
out maximum about 0.4L, and finally con- 
tracted again toward the top. For smaller orifices, the 


spout profiles were similar, except that the constriction 
noted for the in. orifice was absent when the and in. 


187 


Og 
85° CONE 
& Oc =6 
50 Di 


TABLE 


Spout DIAMETER 0.5 AND FT. I.D. 
COLUMNS 


in. 
| | | | 


orifices were used. Thus, the orifice diameter had 
direct effect spout diameter. Data from the 0.5 ft. 
column supported this observation. the variables such 
column diameter, bed height, and air flow which were 
examined, only air flow had easily discernible effect. 
Spout diameter was approximately proportional the 
flow rate above minimum required for spouting, 20% 
increase flow rate causing roughly the same increase 
spout diameter. 


Particle Flow Pattern 


The circulation solids within spouted bed followed 
regular, and easily studied pattern. Particle velocity 
measurements the wall were used estimate solid 
circulation rate. Typical data are given Figure 
velocities measured the circumference 
the column are plotted against fractional bed depth for 
the 0.5 ft. and ft. column. Particle velocity in- 
creased linearly with bed level, and particle velocity 
the wall may represented linear equation havi 
measure the cross flow wheat into the spout. 


Flow Line Correlations 


Measurement particle velocity and direction 
motion the transparent face sectional column 
enables flow lines for particles starting various radii 
constructed. Several flow lines for the ft. column 
obtained this way are shown Figure the top 
half the bed, the flow lines are approximately parabolic, 
the equation which can developed the following 
way. 

Assume the origin the wall, and the top 
the bed. The radius segment then (r,-r) 
distance from the top the bed. Assume also cylindrical 
symmetry flow lines, and uniform bulk density through- 
out the bed. Then material balance 
element volume, with wheat flow the 
top, out the bottom, and from the wall and out 
toward the spout, will result the differential equation 


~ 
u 
0.08 Dc-O5FT 
FT/SEC 
° k = .O57 FT./SEC/FT 
= ° 
0.06 
a 
' 
0.02 
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Figure 6—Wheat flow lines annulus section ft.; 
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Figure 7—Flow line correlation for annular solid flow. 


Integrating with respect the wall 
(5) 


Integrating with respect constant 


definition 


Integration Equations (5) and (7) leads 


and 


where (r,-r,) the radius the spout, and (t-t,) the 
time required for particle starting (r,y,t) reach 
the spout wall Since the particle requires the same 
time fall the distance y,, Equations (8) and (9) 
may equated, and rearranged give 


r(2r, 


where the vertical distance inches from the top 
the bed the point which the particle enters the 
spout, and (r,-r) and are the starting radius and 
spout radius inches, respectively. 


Three sets flow line data are available test this 
given Equation (10). The 0.5 ft. column data 
were obtained photographic and these 
plot linear relationship. The ft. column data 
the higher flow 252 c.f.m. also show linear relation- 
ship about the ft. level, but 210 ie. slightly 
above minimum spouting flow, the 
holds only inches from the top. comparison 


. 


CYCLE TIME 
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SOLIDS 0.2 
CIRCULATING 


WEIGHT PERCENTAGE BED 
Figure 8—Weight distribution vs. cycle time and per cent 


total solids passing each bed level. ft.; 
1.6 ft./see. 


TABLE 


COMPARISON CALCULATED AND OBSERVED SLOPES 
For PARTICLE FLOW ANNULUS 


ft. ft./sec. ft./sec./in. in. in. 


the slopes these lines with ratio shown 
Table The approximate agreement indicates that par- 
ticle velocity data the wall does reflect, least the 
top section the bed, the solid flow pattern within the 
annulus. the bottom half the ft. column, the effect 
the cone flow lines dominant that the simple 
model postulated could not expected apply. 

From the flow line data Figure the time required 
for particles make one cycle can estimated. The 
distribution cycle times for one run shown Figure 
Although some short circuiting wheat occurred 
the top the bed, the major portion the wheat 
90%) took seconds longer make one cycle. 
Also shown the same figure the fraction total solids 
flowing through the spout various levels the bed. 
About two-thirds the wheat entered the spout the 
cone region. Study ft. bed wheat 0.5 
column showed similar solid flow versus height distri- 
bution plot. 


Effect Design Variables Annular Solid Flow 

The effects changing variables such orifice size, 
cone angle, column diameter, bed height, and air flow 
solid flow pattern were evaluated terms particle 
velocity h/L and cross flow Data have been 
summarized Tables and 

The first group data Table show the results 
varying inlet diameter. Orifices from 3/8 in. 1-1/4 in. 
diameter were used the 0.5 ft. column. Total air 
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Vo 
0.6 
x 
4 = 


TABLE 


EFFECT ORIFICE CONE ANGLE 
ANNULAR SOLID FLow 


sec. ft. sec. in. 
1.0 .061 .051 1.4 
| | 


flow was constant. There was significant change the 
total amount solids circulating but there was increase 
the cross flow rate the orifice diameter was increased. 
Thus, more solids entered the core from the cone region 
inlet air velocity was increased. similar series runs 
the ft. column did not show such definite trend. 
Inlet orifices from 1-7/8 in. in. diameter were used, 
but total flow could not held constant because larger 
air flows were required maintain spouting inlet 
diameter was increased. These data show that overall 
circulation rate increased with increasing inlet velocity, 
ie. with decreasing orifice size. 

The second group data show the effect cone 
angle solid circulation pattern. Three cone angles 45°, 
60°, 85° were used 0.5 ft. and ft. column. the 
0.5 ft. column the total flow and cross flow were the same 
for the 60° and 85° cone, and both little less than was 
observed with the 45° cone. the ft. column, the total 
solid flow was about the same for the and 85° cone but 
lower for the 60° cone. Cross flow into the spout appeared 
increase with increasing cone angle. 


Data Table show the effect varying bed height 
and air flow rate solid circulation rate. Since wheat 
does flow into the spout all levels bed, circulation 
rate would expected depend bed height. The 
gas flow required for spouting also increases with bed 
height and these two variables are not independent 
hence the effect bed height not clearly The 
dominating variable affecting solid circulation rate the 
gas flow rate. examination data Table shows 
that increase total air flow results corresponding 
increase particle velocity. 


Air and Solid Flow Correlation 

The relationship between air flow and solid circulation 
rate can developed the following way assumed 
that the pressure drop across spouted bed determined 
the density the wheat-air mixture the spout 
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Figure 9—Relationship between particle velocity and super- 
ficial air velocity. 
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Spout void fraction, can also expressed 


(11) and (12) and rearranging 
APp, 
(13) 


Expressing air and solid flow rates terms superficial 
velocity, putting and taking the area the 
annulus equal the area the column, Equation (13) 


becomes 


This indicates linear relationship 
column velocity and the function given the left 
hand side the equation. Figure which typical data 
for 0.5 ft. and ft. columns plotted, shows that linear 
relationship obtained. Although the ft. column data 
available also show linear relationship, slope 
cept vary considerably from the other column size data. 
Values slope and intercept, the latter expressed terms 
increase slightly with increasing bed depth and 
with decreasing column diameter. Using average values 
for slope and intercept, the following equation relating 
and obtained which fits the 0.5 ft. and ft. 
column data within 15% 


Under minimum spouting conditions 


Lp, Un 


Plotting Lp,/AP results straight line (Figure 
10) with slope 0.14 and intercept 1.0 confirming 
the linear relationship between solid circulation rate and 
air flow and hence the assumption that the pressure drop 
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Figure 10—Relationship between pressure drop and air: 
solid flow ratio minimum spouting. 


across spouted bed determined the wheat-air mix- 
ture the spout. Use Equation (15) (16) along 
with the fact that the pressure drop across the bed 
approximately two-thirds the pressure due the weight 
the bed enables particle velocity and hence total 
solid flow rate calculated for the 0.5 and 1.0 
column. 


Particle Flow the Spout 


The movement particles the spout was clearly 
visible the flat face the semi circular column, and 
was recorded film with mm. Western Fastax cine 
camera using film speed 2000 frames per second. 
These photographs showed that the cone region the 
the cone were constantly collapsing 
particles were being swept the high velocity jet. 
Above the cone region, particles entering the spout appear 
particles with particles the wall the spout. the 
central region the core, however, wheat particles move 
streamline paths with few interparticle collisions. 

Particle velocity different bed levels was estimated 
from the photographs and one set data, minimum 
spouting velocity, shown Figure 11. This based 
particles the centre core, and attempt was 
made determine velocity distribution across the cone. 

The main forces affecting particle motion are: (a) 
accelerating force due frictional drag the particle 
the air stream, (b) decelerating force due gravity, 
(c) decelerating force due collision particles ith 
the wall. This latter force presumably proportional 
the density particles the spout, which varies with bed 
height. first approximation, estimating particle 
velocity, this term will neglected. Using the first two 
forces, the following equation may derived from the 
usual drag Equation (8) 


dt Pp 


Both and are functions the bed height. From the 
air distribution measurements described earlier, can 
estimated, and for the data shown Figure 11, 
ft./sec. the top the bed. The air velocity the 
orifice 80.6 ft./sec. the change spout velocity 
assumed linear with bed height, the relationship for 
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EQUATION 


HEIGHT ABOVE ORIFICE -INS. 


With Equations (17) and (18), and the third relation- 
ship that 


particle velocity various heights can calculated 
successive approximation method. 

The calculated particle velocity curve shown 
Figure 11. The agreement with the observed curve 
fair view the assumptions made. The maximum 
velocity about 10% low, and the curve also predicts 
the maximum somew hat higher the bed. The 
assumption linear relationship Equation (18) 
satisfactory for the upper half the bed, but there 
supporting data for the cone region here the acceler- 
ating force dominant. fact this region, might 
higher than predicted Equation (18) because 
necking down the spout in. was observed. This 
would course increase the magnitude the drag term. 

The velocity the particle leaves the bed can 
estimated from the spout height the equation 


1/2 


The curve based Equation (20) shown Figure 11. 
The particle velocity falls substantially below the observed 
velocity the top ‘of the bed. similar calculation for 
the same bed air flow 252 c.f.m. showed close agreement 
with the observed value. The cause the discrepancy 
210 c.f.m. not known. 

the basis these limited data, particle motion can 
estimated approximately utilizing the frictional drag 
and gravitational forces only, indicating that the deceler- 
ation due interparticle collision small and can 
neglected. 


Conclusions 


The minimum spouting velocity ft. column 


can estimated within 15% accuracy the previ- 
ously proposed Equation (1) which was based 0.5 
ft. columns. 
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suitable general correlation based packed bed 
properties was found for predicting pressure drop across 
spouted bed. Pressure drop was however approximately 
two-thirds the pressure due the weight the bed. 
The pressure drop across spouted bed determined 
the density wheat-air mixture the spout. 

The maximum spoutable bed height can estimated 
from the air velocity through wheat bed which produces 
pressure drop equal the weight the bed. The 
Ergun equation can used for estimating the air velocity, 
and Equation (3) relating minimum spouting velocity and 
bed height can used estimate the bed height. 

The maximum proportion air passing through the 
annulus obtained with large column diameter, large 
orifice size, and small cone angle, and when the spouted 
bed operated near minimum spouting 


The main variables affecting spout size for bed 
wheat are total air flow and column size. Orifice size has 
little direct effect. 

The flow wheat the annulus can estimated from 
particle velocity measurements taken the wall the 
column. Particle velocity data show that about two-thirds 
the wheat enters the spout the cone region, and 
about one-third along the spout. The flow line data within 
the upper half the annulus follows path 
which can derived from simple model. Analysis 
the flow line data also yields data distribution cycle 


times within the bed. 

Column variables such cone angle and orifice size 
affect mainly the rate cross flow solids into the spout. 

The main variable affecting solid turnover rate total 
air flow. relationship between air flow, solid flow and 
pressure drop can derived. 

The principal forces determining particle motion 
the spout are frictional drag and gravitational forces. 
The decelerating force due interparticle collision 
small. The particle velocity can therefore estimated 
from equation based only drag and gravitational 
terms. 


Nomenclature 

Drag Coefficient assumed equal unity 
Constant 
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Level bed, ft. 

Cross flow velocity gradient, ft./sec./ft. 

Bed height above orifice, ft. 

Volume air flow rate, s.c.f.m., 
Radius column, ft. 

Time, secs. 

Linear air flow rate, ft./sec. 

Vertical component particle velocity, ft./sec. 

Horizontal component particle velocity, ft./sec. 
Bed depth measured from top the bed, ft. 
Distance above orifice, ft. 

Density, 

Void fraction 

Viscosity, lb./ft. sec. 

Subscripts 

packed bed 

fluid (air) 

incipient 

inlet 

m => minimum 

wall 

particle 

spout 
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Heats Mixing and 


Equilibrium Calculations' 


BENJAMIN C.-Y. LU? 


the literature, when the available integrated 
forms the Gibbs-Duhem equation are employed for 
extrapolating vapor-liquid equilibrium data iso- 
conditions, the relative partial molal enthalpies 
constant liquid compositions, are assumed either 
investigation the temperature effect heats 
mixing for non-ideal solutions reveals that these as- 
sumptions are invalid and that the quantity varies 
linearly with temperature. Therefore the liquid activ- 
ity coefficients constant compositions should involve 
three temperature terms the form (b/T) 
log for extrapolating purposes. result, more 
general equations are proposed. 


successful application various distillation meth- 
ods chemical industries and the design fractional 
distillation apparatus depend upon the knowledge vapor- 
liquid equilibria the mixture encountered. For this 
reason, most investigations vapor-liquid equilibria are 
made constant pressure conditions. When these data 
are reported, their thermodynamic consistency has usually 
been checked the authors. The purposes testing the 
thermodynamic consistency are indicate that the data 
are self consistent and smooth the data for interpolating 
and extrapolating purposes. One the test methods 
usually employed fit the data means some 
algebraic equations with the condition that the Gibbs- 
Duhem equation the form 


satisfied. The frequently used equations the literature 
are the equations van Laar and Margules rearranged 
Carlson and Colburn the van Laar equation pro- 
posed and Coull Robinson and 
Gilliland and and Coull and other equations 
such presented Wohl Redlich and Kister and 
Benedict, Johnson, Solomon and Rubins 


The applications these equations for extrapolating 
vapor-liquid equilibrium data from one isobaric condition 
another should, however, made with reservations. 


1Manuscript received June 17, 1959. 
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Ottawa, Ont. 
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For example, Yang and van Winkle measured the 
vapor-liquid equilibrium data for the system 
benzene various total pressures. They found that their 


data can represented the van Laar equations, pro- 
posed Robinson and Gilliland 


Xi 2 
Xe 


each the total pressures. The constants and 
Equation (2) are supposed independent temper- 
ature. However, they observed that for the system in- 
vestigated, the values and increase with increasing 
pressure, therefore, the values also increase with increas- 
ing temperature). Eduljee, Kumarkrishnarao and Rao 
correlated the vapor-liquid equilibrium data for the sys- 
tem acetone-water, using form the van Laar equation 
given and Coull and White 


The constants are found functions temperature, 
while they are supposed independent temperature. 
There are many other similar observations the literature, 
indicating that equations such Equations (2) and (3) 
may used for representing the data over the temper- 
ature range certain total pressure, but they are not 
good enough for extrapolating purposes. 


and Coull pointed out that Equations (2) and 
(3) were true, plot vs. 1/T constant composi- 
tion would yield straight line passing through the origin. 
They change activity coefficients with 
temperature for the systems ethanol-water, chloroform- 
ethanol and -water over short temperature 
range and found that when plotted against 1/T 
constant composition, straight line obtained over 
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the range under consideration, but not passing through 
the origin. They obtained 


and suggested the following form the isobaric equation 


for calculating the activity coefficients solution over 
its boiling point range low pressure. Should Equation 
(4) employed extrapolating purpose for the system 
ethanol-water, using the constants obtained and 


Coull the activity coefficients ethanol infinite 
dilution represented Equation (6), 


(5b) 


would always increase with increasing temperature. Otsuki 
and Williams determined the vapor-liquid equilibrium 
data for this system various total pressures from atm. 
300 They observed that the activity coefficients 
ethanol infinite dilution decrease with increasing 
temperature. 


The van Laar equations rearranged Carlson and 
Colburn 


B X2 


represent well the data many binary systems and 
appear the simplest application. Nord (12) sug- 
gested that Equation (7) good perhaps slightly 
better than Equations (3) and (5). known, however, 
that there are many systems which can not represented 
all means Equation (7). This applies even 
binary drocarbon mixtures, whose activity coefficients 
are found functions temperature Since 
Equation (7) does not take into account the variation 
activity coefficients with temperature, any extrapolation 
using the constants obtained one isobaric condition 
another wou'd erroneous most likely. 


Other equations, which involve similar assumptions 
Equations (2), (3), (5) and (7), are usually not ade- 
quate enough for extrapolating isobaric vapor-liquid equi- 
librium data either. The constants obtained one con- 
dition may not used for other conditions. 


seems clear, therefore, that all the integrated forms 
the Gibbs-Durhem equation proposed the literature, 
when they are employed for isobaric equilibrium data, are 
not general enough for extrapolating purposes. The pur- 
poses this investigation are determine the temperature 
effect activity coefficients, using the correlated results 
heats mixing data, and develop more general 
equations for representing and extrapolating binary vapor- 
liquid equilibrium data isobaric 


194 


Temperature Effect Liquid Activity Coefficients 


The activity component, say component 
non-ideal solution defined 
where the chemical potential component Inte- 
grating this equation, 


where refers the molar free energy pure 
component liquid the T,p the solution. Substi- 
tuting the activity coefficient which defined 


into Equation (9), and dividing both sides the resulting 


Partial differentiation this equation with respect 
temperature constant pressure and constant composition 


gives 


= pt R (° (at constant P, 


or 

Rearranging, 


This the expression for the temperature 
the activity coefficients, which molal en- 


thalpy pure component and partial molal 
enthalpy component the solution the same 
temperature. Similar equation may derived for com- 
ponent For simplicity, the dependence activity co- 
efficients pressure neglected the following discus- 
sion. 

The quantity may defined the relative partial 
molal enthalpy component the heat absorbed 
per mole component constant T,p when the pure 
component dissolved large quantity the solution 
without changing the composition the solution. 

seen from Equation (12) that positive quan- 
(heat absorbed) activity coefficient decreases with 
temperature. the other hand, when negative 
quantity (heat evolved), activity coefficient increases 
with temperature. The magnitude variation activity 
coefficients with temperature depends the quantity 
Carlson and inspected the values 
infinite dilution calculated from heats mixing liquid 
the International Critical Tables, and reported values 
high 2000 calories per mole for organic liquid mixtures, 
and 5000 calories per mole for mixtures alcohols and 
water. 


the literature, the integrated forms the Gibbs- 
Duhem equation may separated into two groups. One 
which considers that is, there negligible 
heat effect mixing liquids. Therefore, the activity co- 
efficients are assumed independent temperature. 
The application the integrated forms Equation (1), 
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which are limited isothermal and isobaric conditions, 
isobaric vapor-liquid equilibrium data involves automatic- 
ally this assumption. 


The other group considers constant, constant 
liquid composition. Substituting this expression into Equa- 
tion (12), and integrating, yields Equation (4). 

Equations (2) and (3), the change activity co- 
efficient with temperature constant composition 
considered 


constant 
(In = 
other words, the integration constant Equation 
(4) assumed zero. seen that Equation (13) 
only special case Equation (4), which offers better 


approximation than all the other cases considered. 


Should Equation (4) exact representation the 
variation activity coefficients with temperature, the 
experimental procedure for the determination vapor- 
liquid equilibrium data would very much simplified. 
very unfortunate that there much experimental 
evidence indicating the inadequacy Equation (4). First, 
Equation (4) does not permit the value activity co- 
minimum any temperature, and second, there are indi- 
cations that changes rather rapidly with temperatures. 


Correlation Heats Mixing Data 


The term heat mixing used this paper defined 
the amount heat absorbed the system when 
moles pure component and moles pure com- 
ponent are mixed constant temperature and pressure, 
the amount heat absorbed per mole 
mixture when the pure components are mixed constant 
temperature and pressure, Hence 


A Hy™ = — (Hr): — Me. (Hr pz (14) 


where the enthalpy the solution. Partial differ- 
entiation Equation (14) with respect and gives 


1 


order correlate the relative partial molar enthalpies 
and with temperature and therefore determine 
the variation activ ity coefficient with temperature, the 
knowledge heats mixing essential. The quantity 
the following manner, 


where a(n,,n,), b(n,,n,) and are functions 
and only. Differentiating Equation (16) with 
respect n,, 

ni = 
seen therefore the quantities and have the same 
temperature terms the quantity simplest 


four temperature functions for are discussed here. 


The activity coefficients are independent temperature. 


Therefore constant and 


This the equation proposed and Coull 
Equation (4) represents ‘the condition proposed 
White and Coull and Robinson and Gilli- 
land (16), 


function temperature. 


Therefore 
and 
Case constant composition 
Therefore 
and 


should pointed out, passing, that the constants 
involved Equations (18-23) may either positive 
negative. Therefore, only plus signs are employed these 
equations. 

attempt has been made correlate the quantity 
ture, wtih temperature, using the available data the 
literature. found that most the existing data are 
limited single temperature, making extensive cor- 
relation impossible. However, from the available data for 
non-ideal solutions and some the near ideal solutions, 
the correlated results indicate that the assumptions 


(at constant composition) 
constant (at constant composition) 


involved Equations (1-5) and (7) are too simple 
describe the situation. This explains the inadequacy 
these equations for extrapolating purposes. 


part the investigation may presented here: 


The heat mixing data Bosjnakovic and Grumbt 
(15) for the system ethanol-water (0°-100°C.) are recalcu- 
lated and shown Figure found that low alcohol 
concentration (10-50 wt. alc. 4.2-28.1 mole alc.), 
can represented very well means Equation 
(18), varies linearly with temperature constant 
compositions. However, higher alcohol concentrations, 
the temperature effect apparently more com- 
plicated high temperatures. This may caused 
errors due evaporation condensation, which have 
large heat effects. 


When the relative partial molal enthalpies infinite 
respectiv ely) for the systems: ethanol-water (0°- 100° 
using the calculated Jones, Schoenborn and 
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Figure 1—Heat mixing data for the system ethanol- 
water (15), 


water (0-43.4°C.), carbon disulphide-ethanol (0-15.5°C.), 
ethanol-benzene (0-15°C.) and ether-chloroform (3- 
24°C.), using the calculated values Colburn and Pigford 
(17) are plotted against straight lines are obtained indi- 
cating that they are linear functions temperature. 
part these results shown Figure 


Cheeseman and Ladner investigated the varia- 
tion the heat mixing with temperature for the three 
nearly ideal mixtures: benzene-toluene 
(14.7°-53.3°C.), and observed that the heat mixing 
varies almost linearly with temperature. 

There other evidence the literature supporting 
the applicability Equation (18). 

Scatchard and Ticknor expressed the Scatchard 
excess free energy GF, cal./g-mole, for the systems 
methanol-carbon tetrachloride (25-55°C.) and methanol- 
benzene (25-55°C.) means the following expression 


can readily shown that Equation (24) direct 
support Equation (18). The Scatchard excess free 
energy may expressed 


Since the heat mixing for ideal solution zero, the 
excess enthalpy mixing equivalent the heat 


AH™ 
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Figure 2—Variation partial molal enthalpies infinite 
dilutions with temperature. 


100 


E 


(at constant x)...... (26) 


Differentiating Equation (26) with constant one 
obtains, 


or 


where the integration constant. Therefore 
This agrees with the expression Equation (24). 


and McGlashan measured the heats 
mixing for the systems carbon tetrachloride-cyclohexane 
(10-55°C.) and benzene-diphenyl (10°-60°C.). Their re- 
sults agree with the temperature terms Equation (24) 
and therefore support Equation (18). 


Proposed Equation 


seen from the above discussion that over mode- 
rate range temperature and constant liquid composi- 
tions, the relationship that the quantity linear func- 
tion temperature (Equation (19)) provides much 
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closer approximation than the other relationships pro- 
posed the literature, namely, and constant. 
Therefore, the activity coefficient for component 
constant composition should follow the temperature terms 
indicated Equation (20). The following forms the 
isobaric equation may used for calculating the liquid 
activity coefficient solution over larger temperature 
range than the previously suggested equations the 
literature. 


Redlich and Kister proposed the following equa- 
tion for representing the liquid activity coefficients 
isothermal conditions, 


(31a) 


They have developed these equations expanding the 


function 
temperature. The quantities and Equation (31) 
may modified for representing isobaric data includ- 


ing the temperature terms expressed Equation (20), 


into power series constant 


Therefore 


Equation (33) exact representation Equation (20). 


For systems which can fitted the van Laar 
type equation constant temperature, the following equa- 
tion may employed. 


a+b/T+c log 


a+b/T+clogT 


a+b/T +clogT 


equation somewhat different from Equation (20). 
may seen from actual calculations, however, that 
slight variation the ratio 


does not affect materially the value The advantage 
Equation (34) that only requires three accurately 
determined experimental points (T-x-y). The necessity 
for accurate data, however, cannot over-emphasized. 


and Benedict al. may modified similar man- 
ner introducing Equation (20). 


Although the foregoing discussion limited binary 
the temperature correction terms can readily 
extended into the calculations ternary activity coeffi- 
cients. For example, the equations proposed Redlich 
and Kister for ternary systems constant temperature 
can modified the same manner the modification 
Equation (31), and employed for representing iso- 
baric ternary data. The proposed expression for repre- 
senting liquid activity coefficients function tem- 
perature (Equation 20) may also prove useful the 
liquid-liquid equilibrium calculations. 


Sample Calculation 


The extrapolation the vapor-liquid equilibrium data 
for the system ethanol-water illustrates the calculation. 


The vapor-liquid equilibrium data for the system 


74.79°C. and the heats mixing data for the system 
over the temperature range 0-100°C. are employed for 
determining the temperature effect activity coefficients. 
First all, the thermodynamic consistency the iso- 
thermal data tested means the Redlich-Kister 
relationship (7). 


v2 


x1-0 


observed that the area represented the integral 
Equation (35) equals 0.0081, instead zero. The data 
are then adjusted means procedure similar that 
proposed Broughton and Brearley The observed 
relative volatility corrected using the relationship 


The quantity determined 


log 


The corrected activity coefficients are then calculated. 


isothermal conditions, the van Laar equations 
(Equations (7a) and (7b)) may rearranged give 


eliminating and from these equations. the data 
can fitted means the van Laar equation, straight 
line will result plotting (log vs. (log 
will seen from Figure that the isothermal data can 
not very well represented means the van Laar 
equations. Equation (31) then used fit the activity 
coefficients. The values the constants 


0.5311 
0.1650 
0.0441 


are obtained following the procedure outlined Redlich 
and Kister Hence, 74.79°C., 


and 
log 0.5773(B D/3) 0.298..... (41) 


The heats mixing data Bosjnakovic and Grumbt 
(15) are used obtain the following expressions rela- 
tive partial molal enthalpies, 
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LOG 


Figure 3—Plot (log vs. (log for the system 
ethanol-water 74.79°C, (21), 


0.6 


LOG 


0.2 0.4 0.6 
activity coefficients and experimental 


values for the system ethanol-water 25°C. 
experimental (2°) 


859.8 2.950 T........... (44) 
1238.7 3.400 T........... (45) 


should mentioned that these expressions are obtained 
giving emphasis low temperature values. 


Substituting Equation (42) into Equation (12) and 
integrating 


198 


Figure activity coefficients and experimental 
values for the system ethanol-water 
Calculated 
experimental (11) 
experimental 


19.217 log 
the integration constant determined from Equation (39). 
Similarly, 
and 


The quantities and expressed Equation (32), 
may determined from these expressions. Therefore, 


and 


Equation (49) used calculate the liquid activity 
coefficients for the system ethanol-water atm. and 
25°C. The calculated values are compared with the ex- 
perimental data Figures and both cases, good 
agreement obtained. 
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(22) 
(23) 


The removal droplets from dioctyl phthalate 
mist fluidized bed silica gel was studied. The 
mist used had narrow range droplet diameters, 
which averaged 0.87 microns size. Silica gel par- 
ticles were the range —150 mesh +200 mesh. 


was found that removal efficiencies were only 
slightly affected 50% change inlet mist con- 
centration, 300% variation depth the 
fluidized bed. Collection efficiencies varied from 
about 90% low superficial gas velocities 70% 
high velocities. 

shown from these results that collection 
primarily particle diffusion. Results can cor- 
related terms equations for mass transfer 
fluidized beds. The possible uses fluidized bed 
collecting device are discussed. 


use fluidized bed remove droplets from mists 
been reported only once the literature 
Meissner and Mickley this work used such bed 
remove droplets the size range 2-10 microns diameter 
from sulphuric acid mist. Over the range conditions 
used, was found that the removal efficiency was inde- 
pendent mist concentration, and increased somewhat 
bed weight increased. Efficiencies increased with in- 
creased superficial velocity through the bed. Preliminary 
work with ammonium nitrate dust was also reported, and 
the above results were confirmed with the exception the 
effect velocity the degree removal. 


From the results obtained Meissner and Mickley, 
the fluidized bed seems offer number advantages 
aerosol collection device. would have nearly 
constant pressure drop, independent the age the bed. 
Collection efficiencies can well over 90%, and these 
have also been found independent bed 
ations inlet mist concentrations seem have effect 
collection performance over wide ranges, 
ations the level solids the bed are minor 
importance also. 

Fluidized beds porous solids (for example, spent 
cracking catalyst, activated alumina, and silica gel) have 
own weight liquid before bed agglomeration due 
saturation causes any effect collection efficiency. This 
was also true the collection ammonium nitrate dust. 


1Manuscript received May 30, 1959. 
2Associate Professor, Dept. Chemical Engineering, The University 
British Columbia, Vancouver, B.C. 

8Hooker Electrochemical Limited, North Vancouver, B.C. 

Contribution from the Department of Chemical Engineering, The Uni- 
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Such bed has limitations with respect temper- 
ature. could operated continuously, and the bed 
solids regenerated and the aerosol material recovered 
destroyed auxiliary unit. 

The most serious limitation would seem that this 
device primarily suitable for mists, although appears 
that could expected perform equally well for 
some dusts. Pressure drops through the bed are about 6-8 
inches water, which greater than that through 
other devices comparable efficiency. 

the present work, the performance fluidized bed 
silica gel particles was investigated for the removal 
homogeneous dioctyl phthalate mist. The average drop- 
let diameter all runs was 0.87 microns, with all particles 
the size range 0.5 1.1 microns. 

sensitive colorimetric method was developed for the 
analysis dioctyl phthalate. This was necessary because 
was found that traces silica gel dust samples were 
sufficient cause erratic results when analysis ultra- 
violet absorption was attempted. 


Theory 


The mechanism collection fluidized bed could 
easily visualized occurring either impingement due 
inertial effects, droplet diffusion, combination 
both. The ranges conditions which these mechan- 
isms give effective contribution the collection effici- 
ency have been investigated Albrecht Langmuir 
Ramskill and Anderson and Gillespie (5), connection 
with the theory aerosol filtration. has been shown 
Ramskill and Anderson that there certain particle 
size for each set filtration conditions which will result 
minimum collection efficiency; this efficiency increases 
particle size either increased decreased. For many 
combinations filter characteristics this “most difficult” 
size lies the range 0.2-2.0 microns diameter. 
shows further that correlation collection efficiencies with 
filter characteristics made much more difficult the 
fact that many interceptions aerosol particles col- 
lecting surfaces are non-effective, particularly high 
velocities. 

Meissner and Mickley believed the major collection 
mechanism their work was due collision effects, 
although the small effect bed depth could not ade- 
quately explained that basis. Their results were repre- 
sented equation the form 


(1) 


where varied from 0.16 0.34. This equation indicates 
the dependence collection efficiency gas velocity, 
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and the number collecting particles. sulphuric 
acid mist was used this work, having average droplet 
diameter from 2-6 microns. The collecting solids were 
various size distributions. 


The inertial effect has been shown proportional 
Dm? Vo/Dp, alternatively, the number 
ency” estimate this value for the results 
Meissner and Mickley indicates that for all the bed solids 
used the inertial effect large, and therefore should 
constitute the primary collection mechanism. The in- 
crease collection efficiency with increased gas velocity 
and decreased size collecting particle, which was ob- 
served this work, also confirms this assumption. 

collection occurred diffusion mechanism, the 
simplest approach would repiesent the process 
equation the same type those applied mass 
transfer operations. For uniformly sized droplets, the 
decrease mass concentration would proportional 
the number droplets collected. The driving force for 
droplet diffusion would the mass concentration, inas- 
much the process not reversible. The mass transfer 
equation could then written 


where the collecting surface per unit volume bed 
solids, and “mass transfer” coefficient for particle 
diffusion. The incremental bed depth can expressed 
where the bulk density the bed, and 
its cross sectional area. Further, can expressed 
Am, the collecting surface per unit weight solids. 
Making these substitutions: 


CdW. 


Integrating 


The dependence the mass transfer coefficient and 
effective collecting surface bed depth, well 
variations the degree fluidization are difficult 
estimate. the usual assumption made that the mass 
transfer coefficient depends primarily the flow rate for 
given system then the above equation can integrated. 
interesting note that diffusional process the 
principal collecting mechanism, the collection efficiency 
should decrease the flow rate through the bed in- 
creased. This feature has been discussed Ramskill and 
Anderson (4), 


Experimental Apparatus 


The mist used was generated Lamer-Sinclair 
homogeneous aerosol generator, shown Figure using 
phthalate the liquid. portion dry filtered 
air was metered and passed through nucleation chamber, 
where nuclei were supplied from heated coil freshly 
coated with pure sodium chloride. The bulk the air 
was metered, passed through preheater, and into the 
vaporizer, containing the dioctyl phthalate. The mixed 
vapor and nuclei entered the reheater, and then were 
slowly cooled the double walled chimney, 


portion the mist laden air stream could bled 


off dilution air could added before the 
mist entered the fluidized bed, The fluidized bed used 
was constructed 2-inch inside diameter glass tube 
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fitted with outlet cyclone and collector for dust re- 
moval. Bed solids used were —150 200 mesh silica gel 
supported 200 mesh screen. 

Sampling the inlet outlet mist was done passing 
the entire flow for suitable period time through 
filter holder and wet test meter. The filter used was 
Millipore aerosol filter. 

negligible amount mist was collected the upper 
empty part the column and the cyclone. One filter 
was found adequate collect all the mist sample 
air. 

Particle size the aerosol generated was determined 
periodically collection samples coated slides, and 
microscopic examination. 


Analaysis Dioctyl Phthalate. 


The analytical method was based the reaction be- 
tween organic esters and hydroxylamine alkaline solu- 
tions form droxamic acid. The hydroxamic acids 
give colored complexes with metals, notably with iron 
compounds. This method was first quanti- 
tative method for esters fatty acids, and 
was applied the present work aromatic esters. Sub- 
sequent this work more detailed investigation the 
reaction was published Goddu, and Wright (). 

aerosol sample giving 0.1 1.5 mg. dioctyl 
phthalate was collected Millipore filter. This was 
extracted with ml. petroleum ether, the solution 
was added 0.3 ml. 2.5% solution NaOH 95% 
ethanol. The mixture was evaporated dryness 
water bath 67°C., removed and cooled 25°C., and 
ml. solution ferric perchlorate and perchloric 
acid ethanol added. After standing for few minutes, 
the transmittancy was compared blank 540 mu. 

Amounts dioctyl phthalate could determined 
with accuracy better than +0.02 mg. this method. 
Typical calibration charts for three organic esters are 
shown Figure 


Results and Discussion 


The mist used all runs had average particle 
diameter 0.87 microns, with about 95% all particles 
lying the range 0.60 1.1 microns. Experimental 
efficiencies were measured averaging several deter- 
minations both inlet and outlet concentrations. 

Figure the effect inlet mist concentration 
removal efficiency shown. Although the range con- 
centrations not wide, apparent that 50% change 
inlet mist loading has effect the collection effici- 
ency observed. This agrees with the results reported 
previously 

The variation collection efficiency with superficial 
velocity shown two typical plots Figures and 
the bed depth used one case being double that the 
other. Examination these data shows that the effect 
bed weight mist removal small, and this work was 
less than the experimental error determining efficien- 
cies. Therefore, the dependence the efficiency, C,/C,, 
bed depth indicated either Equations (1) (2) does 
not appear present, again agreement with previous 
work Several tests varying velocities were carried 
out with the inlet screen place but with the bed empty. 
About 15% 20% the mist was collected passage 
through the inlet section and screen. However, these 
are integral parts the fluidized bed, attempt was 
made correct efficiencies this amount. 

can seen from Figures and that the removal 
efficiency decreases the superficial gas velocity 
the bed increases. The opposite effect was found 
Meissner and Mickley. collection was primarily 
collision mechanism due inertial forces, would 
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CONSTANT TEMPERATURE BATHS 


Figure 1—Sketch the apparatus. nucleation coil. 

preheater. vaporizer. reheater. cooling chimney. 

bleed-off line. dilution line. fluidized bed. mist 
samples. gas meter. 


expected that efficiencies should improve with increase 
the gas velocity through the bed. However, collec- 
tion diffusional mechanism due Brownian motion, 
then the observed effect velocity expected. 


estimate the relative importance inertial col- 
lection can found calculating the “separation num- 
ber”, for the droplets and bed particles used, 
and determining the “target efficiency” which might 
expected for single bed particle use the chart 
prepared Langmuir and Blodgett The value the 
separation number for the work presented here can 
expressed Vo. apparent that for the range 
values given that the separation number well 
below the value about 0.04, the theoretical limit below 
which collection impingement single sphere 
would zero. Even allowing for the dissimilarity 
behavior between single spheres and fluidized bed 
spheres, indicated that collection inertial methods 
cannot the major mechanism, and unlikely con- 
tribute even significant fraction the total collected. 
comparison the magnitude the separation numbers 
shows that inertial effects the work Meissner and 
Mickley were the order 200-500 times great 
these experiments. 

considering diffusional mechanism, expressed 
Equation (2), the dependence the mass transfer 
coefficient, kp, the gas mass velocity can estimated 
from correlations for fluidized beds. Use the equations 
proposed Chu, Kalil and Wetteroth would indicate 
that should proportional for the range 
Reynolds numbers covered these tests. 


This would imply that the ratio, should vary 
The results given Figure are plotted 
this way Figure show the degree correlation 
obtained. While the resulting plot satisfactory, the data 
are not sufficiently extensive confirm value 0.78 
the most suitable exponent. The point where fluidization 
developed shown very clearly plot this type. 
Further results have indicated similar transition high 
superficial velocities when “slugging” occurs the bed. 


Venturi Collector—As matter interest, and check 
the ability predict collection mechanism use the 
relations for inertial effect, another type device was 
used which inertial effects would expected 
much more significant. This consisted small glass 
venturi tube set over the inlet pipe the fluidized bed 
that all the gas flow was through the venturi, which 
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Figure Analysis mist concentrations. dioctyl 
phthalate. n-butyl phthalate. ethyl benzoate. 


was turn surrounded the bed silica gel 
Two holes were made the throat the venturi through 
which bed solids were aspirated. The gas-solid mixture 
left the venturi above the solids level that the solids fell 
back the bed. such device, high relative velocity 
between mist droplets and solid collecting particles can 
obtained the throat the venturi. 

Throat velocities were varied from about 300 cm./sec. 
3000 cm./sec. Solids aspiration was such that about 
mist droplets were present for each solid collecting 
particle, except the highest velocities. Over the above 
range velocities the collection efficiency increased 
slightly with throat velocity, from about 37% 42%. 
According Langmuir collection single sphere 
over this range should vary from about 10% 58%. 
view the assumptions the theory, and its application 
this case, the agreement surprisingly good, and 
illustrates the point that one can predict qualitatively the 
nature the collection mechanism 
device this means. 
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Figure 3—Effect inlet mist concentration. 
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Figure 4—Variation removal efficiency with superficial 
gas velocity through bed. 


Figure 6—Correlation removal efficiencies and superficial 
velocities, bed wt. gm./sq.cm. fluidized bed. 
packed bed. 


Summary 


fluidized bed fairly efficient device for removing 
mist droplets size range which usually presents serious 
collection difficulties, that from 0.5 1.5 microns 
diameter. The collection mechanism appears prim- 
arily due diffusional effects. Collection efficiencies are 
not affected seriously variations inlet mist loadings, 
bed age, the amount bed solids, within fairly wide 
limits. 


possible predict qualitatively from theory, col- 
lection behavior fluidized bed gas flow increases. 
Additional data would make this prediction quantitative. 
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Velocity, cms./sec. 


Figure 5—Variation removal efficiency with superficial 
gas velocity through bed. 


very possible that fluidized beds could have 
applications mist dust collecting devices, particularly 
for those droplet sizes usually considered form the most 
stable aerosols. 


Nomenclature 

collecting surface per unit weight fluidized solids. 
collecting surface per unit volume fluidized solids. 
Mist concentration, weight per unit volume. 

Diameter Mist droplet. 

Diameter the collecting solid. 

gravitational constant. 

constants. 

cross sectional area fluidized bed. 

terminal velocity mist droplet. 

superficial gas velocity a.fluidized bed. 

weight bed solids per unit cross sectional area bed. 
depth fluid bed. 

bulk density bed solids. 

density mist droplet. 
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Fluidization Beds Rice and Wheat 


LIANG-TSENG FAN and CHARLES SWARTZ 
Engineering Experiment Station, Kansas State University, Manhattan, Kansas. 


fluid flowing upward admitted into the 
bottom bed granular material, pressure 
drop through the bed will result. This pressure drop will 
increase with increase fluid rate until point equili- 
brium reached which the pressure drop through the 
bed multiplied the cross-sectional area the column 
equal the weight the bed and which the fluid- 
like motion the bed commences The fluid velocity 
corresponding this point defined onset minimum 
fluidization velocity the bed. 


new generalized correlation, which relates such 
factors fluid viscosity, fluid density, particle density, 
shape factor particle and void fraction onset point 
the onset velocity fluidization, was recently proposed 
Leva The data used formulating this cor- 
relation were mostly obtained from experiments using 
metallic mineral solid particles with small diameters. 
The validity the correlation fluidization cereal 
grains, which mostly consist relatively coarse granules 
with characteristics quite different from metallic and 
mineral particles, has never been tested. is, however, 
known that when heat being transferred from the 
column wall into the bed, fluidization onset accom- 
panied sharp increase heat transfer coefficients 
Klassen and Gishler recently published the results 
their fluidization heat transfer study which includes such 
cereals rice and wheat. Since their data include runs 
made during the packed range, the plot heat transfer 
coefficients against the flow rate should show breaking 


point the onset fluidization. Application the 
generalized correlation Leva, the fluidization 
conditions Klassen and Gishler yields the fluidization 
onset velocities 490 and 700 for 
rice and wheat respectively. 

adapting the generalized correlation beds 
particles with irregular shape, necessary have some 
means determining average diameter the 
The recommended procedure for finding the diameter 
narrow cut solids take the geometrical average 
adjacent sieve openings and However, the sizes 
relatively large grains such rice and wheat are gener- 
ally defined and measured terms three dimensions, i.e. 
length, width and thickness Comparison the results 
calculations using the square root the product 


two smallest dimensions width thickness) gave results 
which are mentioned above and are closest the experi- 
mental results approximately 500 for rice 
and approximately 700 for wheat. 

further verify the generality the finding, several 
sets experiments were conducted using pyrex glass pipes 
various sizes the fluidizing columns. The conditions 
and results the experiments are summarized Table 
and Table respectively. cases did the deviation 
the calculated values—using the geometrical mean width 
and thickness the grains the mean diameter the 
grains—from the experimental data exceed 5%. 

The results the present work clearly indicate that 
the generalized correlation proposed Leva, for 


TABLE 
EXPERIMENTS 


Run Temp. Press. Particle Particle Size Column 

No. atm. Fluid Particle Density in. in. Diameter 
in. 

70.0 1.0 air rice 91.4 0.205 0.0900 0.0660 
76.0 1.8 air wheat 77.4 0.241 0.100 0.105 
73.1 air wheat 77.4 0.241 0.100 0.105 
70.0 air wheat 74.3 0.229 0.0958 0.0924 
70.0 air wheat 74.3 0.229 0.0958 0.0924 
TABLE 


EXPERIMENTS 
ONSET VELOCITY 


Run Calculated Observed Deviation 
No. 
646 +3.03 
900 870 —3.33 
935 920 
930 890 —4.30 
925 900 —2.70 


the onset velocity fluidization may applied fluidi- 
zation rice and wheat gas when the proper dimen- 
sions the grain (square root product two dimen- 
sions, width and thickness) are used estimating the 
mean diameter the particle. 
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techniques and growing consumer 
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HEN fluid flowing upward admitted into the 

bottom bed granular material, pressure 
drop through the bed will result. This pressure drop will 
increase with increase fluid rate until point equili- 
brium reached which the pressure drop through the 
bed multiplied the cross-sectional area the column 
equal the weight the bed and which the fluid- 
like motion the bed commences The fluid velocity 
corresponding this point defined onset minimum 
fluidization velocity the bed. 

generalized correlation, which relates such 
factors fluid viscosity, fluid density, particle density, 
shape factor and void onset point 
the onset velocity fluidization, was recently proposed 
Leva The data used formulating this cor- 
relation were mostly obtained from experiments using 
metallic mineral solid particles with small diameters. 
The validity the correlation fluidization cereal 
grains, hich mostly consist relatively coarse granules 
with characteristics quite different from metallic and 
mineral particles, has never been tested. is, however, 
known that when heat being transferred from the 
column wall into the bed, fluidization onset accom- 
panied sharp increase heat transfer coefficients 
Klassen and Gishler recently published the results 
their fluidization heat transfer study which includes such 
cereals rice and wheat. Since their data include runs 
made during the packed range, the plot heat transfer 
coefficients against the flow rate should show breaking 


point the onset fluidization. Application the 
generalized correlation Leva, the fluidization 
conditions Klassen and Gishler yields the fluidization 
onset velocities 490 and 700 for 
rice and wheat respectively. 

adapting the generalized correlation beds 
particles with irregular shape, necessary have some 
means determining average diameter the particles. 
The recommended procedure for finding the diameter 
narrow cut solids take the geometrical average 
adjacent sieve openings and However, the sizes 
relatively large grains such rice and wheat are gener- 
ally defined and measured terms three dimensions, 
length, width and thickness Comparison the results 
calculations using the square root the product 
two smallest dimensions width thickness) gave results 
which are mentioned above and are closest the experi- 
mental results approximately 500 for rice 
and approximately 700 for wheat. 

further verify the generality the several 
sets experiments were conducted using pyrex glass pipes 
various sizes the fluidizing columns. The conditions 
and results the experiments are summarized Table 
and Table respectively. cases did the deviation 
the calculated values—using the geometrical mean width 
and thickness the grains the mean diameter the 
grains—from the experimental data exceed 

The results the present work clearly that 
the generalized correlation proposed Leva, for 


ABLE 


CONDITIONS OF EXPERIMENTS 


Run Temp. Press Particle Particle Size Column 
No. atm Fluid Particle in. in. Diameter 
76.0 1.8 air wheat 0.241 0.100 0.105 
2.2 air wheat 77.4 0.241 0.100 0.105 
70.0 air wheat 74.3 0.229 0.0958 0.0924 
70.0 air wheat 74.3 0.229 0.0958 0.0924 
TABLE the onset velocity fluidization may applied fluidi- 


RESULTS EXPERIMENTS 
ONSET VELOCITY 


Calculated Observed Deviation 


Ib. /hr.-ft.2 Ib. 
900 870 
935 920 
930 890 
925 900 


zation rice and wheat gas hen the proper dimen- 

sions the grain (square root product two dimen- 
sions, i.e. Width and thickness) are used estimating the 
mean diameter the particle. 
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